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ABSTRACT 

Often a pattern recognition system is broken into two parts “ feature 
selection and classification. This note deals with feature selection. 

In 1971, the Earth Resources Laboratory (ERL) received a version of the 
Purdue pattern recognition software (LARSYSAA) . ERL was preparing to pro- 
cess the very large quantities of data as output by the Bendix 24 channel 
scanner, Eppler*s Digital Table Look— up Classifier (DTL) was being im- 
plemented at ERL primarily because of its speed advantage . The Purdue 
feature selection program was quite useful but was limited to twelve 
features (or channels) and did not select channel subsets in a manner sug- 
gested by Eppler to be used with his table look-up classifier. He sug- 
gested using a different subset of channels for each material (class), i.e., 
the subset which is best suited for detecting a particular material. 

Feature selection software has been developed at ERL that is capable 
of inputting up to 36 channels and selecting channel subsets according to 
several criteria based on divergence. One of the criterion used is compatible 
with the table look-up classifier requirements. The software indicates which 
channel subset best "separates" (based on average divergence) each class from 
all other classes. The software employs an exhaustive search technique, and 
designing the software such that computer time did not become prohibitive 
was a major goal. A typical task to select the "best" 4 of 22 channels for 
12 classes takes 9 minutes on a Univac 1108 computer. 



Table of Contents 


Contents 

Introduction 1 

Feature Selection 2 

Divergence 2 

Criteria 5 

Computation Time 7 

Core ^ 

Uses 

Summary H 

References 13 

Appendix 15 

Example of Output 16 

Software Description 30 

Functional Description 31 

Logic Flow 32 

Univac 1108 Exec 8 Run Stream 33 

Fortran Listings 34 

Example of Input Data 64 

Description of Input Data 


69 



INTRODUCTION 


A pattern recognition system is often broken into two parts - feature 
selection and classification. This note deals with feature selection. 

(For a much more detailed discussion of pattern recognition and feature 
selection, see [1, 11, 12].) In 1971, the Earth Resources Laboratory 
(ERL) received a version of the pattern recognition software LARSYSAA. 

This software was originally developed by the Laboratory for Application 
of Remote Sensing (LARS) at Purdue University. The system is well-known 
and is described by Landgrebe [8] and Ratcliff [16]. 

ERL was preparing to process the very large quantities of data as out- 
put by the Bendix 24 channel scanner [19]. Eppler's Digital Table Look- 
up Classifier (DTL) was being in^lemented at ERL primarily because of its 
speed advantage [4]. The LARSYSAA feature selection program ($SELECT) was 
quite useful, but it was limited to 12 channels and did not rank the channel 
subsets in a manner suggested by Eppler, i.e., by the four channels which 
were best suited for detecting a particular material (class). It was 
necessary to develop feature selection software which would handle more than 
12 channels and would rank them in a manner more compatible with the DTL. 

This software- CHOICE - is coded completely in Fortran V and runs on a 
Uni vac 1108 conputer. The maximum number of classes and channels are cur- 
rently set at 32 and 24 respectively and the program uses less than 32K, 
i.e., 32,000 Univac 1108 computer words. These limits are not fixed. 

More details can be found in the discussions under Computation Time and Core. 

/ 
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FEATURE SELECTION 


A brief description of the notions of pattern recognition and feature 
selection as they are often applied in processing multispectral scanner 
data follows . 

A series of n measurements are made on an object. This series could 
be composed of measurements of intensity of spectral radiation in n bands 
(channels) of a multispectral scanner. A series of measurements on an ob- 
ject can be considered an n-dimensional observation vector. Given m multi- 
variate normal populations (classes) , we are to decide to which class the 
observation vector most likely belongs. In general* it is not practical 
to use all n measurements in making this decision (classification) because 
of computer time constraints. We want to select some subset k of n channels 
which will enable us to classify the data in an accurate manner. Often 
there are many subsets to choose from. For example, the DTL requires that 
k=4 and in working with 24 channel scanner data, there are more than 10,000 


possible channel subsets 
channel subset . 



We use divergence to help us choose a 


DIVERGENCE 

All of the criteria used by CHOICE are based on divergence. The appli- 
cation of divergence to feature selection was proposed by Marill and Green 
[10], and an interesting discussion on the subject is given by Fu [11]. 
Divergence is given a very general definition by Kullback [?] and is based 
on considerations from information theory. Kullback argues that divergence 
is an appropriate "distance” between arbitrary distributions and asserts it 
to be "a measurement of the difficulty of discriminating" between two dis- 
tributions [7, 10], In the case of two multivariate normal populations with 
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equal covariance matrices, it can be shown that the probability of classi- 
fication error is a monotonically decreasing function of divergence [11] . 

In the less restricted case of more than two populations with unequal co- 
variance matrices, there is experimental support for using divergence [10] . 

Criteria based on divergence have come under attack because they have 
not been explicitly expressed in terms of probability of error in the more 
general cases. Fuirther, specific examples can be cited where the average 
pairwise divergence criterion does not yield the minimum probability of 
error [6, 15]. Although these criticisms are valid, criteria based on 
divergence have proved to be a useful guide in working with multispectral 
scanner data and in attempting to improve classification results. 

In the case of several multivariate normal populations (classes) , 
the pairwise divergence between class i and class j reduces to: 

= h trj(K^- Kj) h tr (u^- uj)^(Kj^+Kj^) (u^-Uj) (1) 

where : 

tr - trace 
T - transpose 

ui - kxl mean vector of class i 
Ki - kxk covariance matrix of class i 

- inverse of the covariance matrix of class i 

CHOICE actually computes the quantity j , where . 

Figure 1 illustrates some simple exati 5 )les of class pairs and their 
associated divergences. The data from which these statistics were calcu- 
lated were acquired by an aircraft borne 12 channel multispectral scanner. 
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Figure 1. Statistics calculated from actual data (flight line C-1)* Note the relation 
between a large value of pairwise divergence and good separability* 





This flight line (Flight line C-1) was flown June 28, 1966 for LARS and has 
been widely reported on [2, 9]. In Figure 1, probability density functions 
[fi(x)l are plotted for four classes in the single channel case (channel 1). 
The mean response of channel 1 for soybeans (170,0) is the average reading 
of this channel as it "looked at" those soybean fields which were used in 
training the pattern recognition algorithm. The channel readings for a 
class are usually assumed to be normally distributed about the class s 
mean response. Thus, given the mean response and the standard deviation of 
these responses for soybean training fields, the probability density function 
for soybeans can be graphed. Figure 1 illustrates, and the pairwise diver- 
gence (YA = 3.1) indicates, that channel 1 would not be a "good” channel 
to discriminate rye from alfalfa. However, it would be a "good" channel to 
separate soybeans from red clover as evidenced by the figure and by SR=56.9. 

CRITERIA 

Assume that we are given n channels and are asked to choose a subset 
k of these, where k<n, by some criterion. CHOICE ranks the channel sub- 
sets by five separate criteria* The first two of these are the same as 
those available with $SELECT. Unfortunately, there is no one criterion that 
is always "best". Thus, the software presents several rankings to the in- 
vestigator, An example of these channel subset rankings with k=4 for each 
criterion is given in the Appendix under Example of Output . Data from 
Flight Line C-1 was used in producing this output and is also listed in the 
Appendix under Example of Input Data. Following is a description of the 
fiveT criteria used by CHOICE: 

Criterion 1: Maximizing the arithmetic mean of the pairwise divergences. 

The average pairwise divergence for a given channel set can be computed by 
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-m 


Div 


m(m-l) j=i i=j+l 


X) ^ > where m is the number of classes 


This average is computed for all possible subsets. The channel 
subsets which yield the largest averages are saved and printed. 


Criterion 2. Maximizing the minimum pairwise divergence. 

The minimum pairwise divergence for a given channel set is given by 

Div^in = min {Otj ( . j=i+l* m 

The minimum is computed for all subsets. Those subsets which produce 
the largest ininimums are saved and printed. 

Criterion 3. Maximizing the average divergence by class. 

This is the criterion used at ERL most often in conjunction with the table 
look-up. Since the DTL can be run with a different set of channels for 
each class, this criterion is used to help decide which subset of channels 
best separates a class from all other classes. The average divergence for 
class j for a given channel set is computed by 

“^''avg (class 3) ° 

The class average is computed for all classes for each channel set. This 
j[^g repeated for all possible channel sets • The channel subset which yields 
the largest average for each class is saved and printed. 

Criterion 4. Maximizing the geometric mean of the pairwise divergences. 

It seems reasonable to expect the difference between a divergence of 10 
and 30 to be more important than the difference between a divergence of 50 
and 70 even though the delta in both cases is 20. For example, in the one- 
dimensional two-class case with variances equal to one, we get the following: 
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10 

30 

50 

70 


Average Probability of Classification Error 


13.14% 

2.62. 

.62 

.15 


Criterion 4 gives more weight for an improvement in separation from 10 to 
30 than for an improvement from 50 to 70. The more widely available 
criterion 1 considers these improvements of equal value. The geometric 
mean is computed for each channel subset as follows i 
1 m-1 m 

“i^'geo = TT “iJ 

J=1 X=J+1 

This mean is con^juted for 
subsets which yield the largest geometric means are saved and printed. 

Criterion 5. Maximizing the pairwise divergence for each class pair. 

DjQj is computed for each subset of channels . That channel set which 
yields the maximum divergence for each class pair is saved and printed. 

This is most useful in working with hard— to— separate classes. In the 
worst case, it will indicate to the investigator that no matter which 
channel set is selected, it will not be possible to separate the two classes. 
In less severe cases, it will indicate the most separability to be expected 
and which channel set will provide it. (This assumes that the channels will 
not be transformed. See discussion under Uses and [14],) 

COMPUTATION TIME 

There has been some discussion that it is not feasible to do an ex- 
haustive search when the number of channels is large [2, 17]. By incor- 
porating some efficiences, we have found we can routinely process 24 


all possible subsets 


- 0 - 


and those channel 
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channel scanner data using the exhaustive search technique. A typical task 
to select the "best" 4 of 22 channels for 12 classes takes 9 minutes on a 
Univac 1108 computer. 

In order to handle more than 12 channels in a timely manner, we 
decided it would be necessary to rewrite the feature selection software and 
not just modify LARSYSAA. There were two key changes which accounted for 
considerable time-savings. 

1. Input /Output 

No input/output (I/O) is done during the exhaustive search loop. 

Because of their generality, Fortran formated read and write commands are 
quite expensive in both throughput and in execution time. CHOICE circum- 
vents doing an I/O after each set of divergence calculations by immediately 
testing whether this channel set ranks in the top group of channel sets for 
each criterion. Each group is made up of those sets which* have exceeded 
all the preceeding channel sets in maximizing a particular criterion, if 
the set is not ranked in any top group, it is discarded. If it is top- 
ranked, then the channel set is saved using a bit manipulation function. 
With this technique only one 36 bit Univac 1108 word is required to store 
a channel set. (Hence the limitation of 36 channels.) At this point the 
pairwise divergences are discarded in order to conserve core. The search 
is continued until all possible subsets have been considered. At the com- 
pletion of the search, the divergences for those top-ranked channel sets 
(which usually number less than 50) are recomputed and printed. 

2, Matrix Inversion 

In a case where we have 12 classes and are asked to select the best 4 
of 22 channels, we must perform over 85,000 matrix inversions, and these 
comprise a significant portion of the total computation time. A procedure 
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based on the bordering technique [5] was incorporated into CHOICE to invert 
the lower triangular matrices. This procedure takes advantage of the fact 
that the matrices are symmetric. In a test, this procedure was about five 
times as fast as the more general and widely used matrix routine MINV. 

MINV uses the Gauss-Jordan reduction technique and does not capitalize on 
the symmetry of the matrices. 

An empirical formula which can be used to roughly estimate the computa- 
tion (CPU) time of CHOICE in minutes is: 

time = 

50,000 

k - number of channels in subset to be chosen 
m - number of classes 

Take, for example, the theoretically important two-class case. Let us 
estimate the time to choose the best 4 of 36 measurements. 

50,000 

time = 1.2 minutes 



where n - total number of channels 
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Parameter 


Meaning 


Current Setting 


s 

Pi 


P2 


maximum number of show requests 


10 


maximum number of top-ranked sets 

to print for criteria 1 and 2 20 

maximum number of top-ranked sets 

to print for criteria 3 10 


The following formula can be used to calculate the amount of core 
required (In U1108 words) . 

Core = 10,500+m[n+n(n+l)/2+k+k(k+l)/2i4]+k(s+l)+n(n-l) (px+4)/2+4p2 

All parameters must be greater than or equal to one* The number of channels 
oiust be less than or equal to 36* and the number of classes must be greater 
than or equal to 2. 


USES 


Just as with the LARSYSAA separability processor, there are several 
areas in which the investigator may find CHOICE useful in processing multi- 
spectral data. 

1) The primary purpose of the program is to indicate which subset of 
measurements might best be used in deciding to which of several 
classes an unknown element belongs . 

2) The program can indicate those classes which cannot be separated 
with the given set of measurements, no matter which subset we choose. 
For example, CHOICE may indicate that white sand and concrete look 
spectrally very similar when observed from space under certain illu- 
mination conditions. The investigator must decide whether to group 


show request refers to a request by an Investigator to see the pairwise 
divergence printed for a channel set of interest — one that is not necessarily 
top-ranked. 
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these materials into one class or to drop the materials from the 


classification . 

3) Training fields for the same class can often be quite different 
spectrally. A decision must be made as to which training fields 
can be pooled into one population and which cannot. CHOICE can aid 
in this decision by indicating whether the two training fields are 
spectrally "close" or not. "Closeness" would be indicated by 
relatively low pairwise divergences . 

4) Recently there has been developed a feature selection technique 
that seems very promising [13, 14]. Basically this technique 
determines a linear transformation Sx which can be used to reduce 
the dimension of the data from n to k where n is greater than k. 

This technique determines a k by n matrix B which maximizes the 

so called B-average divergence. (One disadvantage of this tech- 
nique is that it is currently limited to 9 classes.) It is re- 
commended that the best k of n channels be chosen as the initial 
guess for the B matrix [3] . This is done in order to increase the 
probability that the maximum found iteratively will be the global 
niaximum. It seems that CHOICE could be quite useful in initializing 
B, especially when n is greater than 12. 

SUMMARY 

In this note, feature selection software (CHOICE) was described. CHOICE 
was developed mainly to enable processing 24 channel scanner data and to 
rank the channel subsets in a manner more compatible with the table look-up 


11 



classifier. A brief justification of the use of divergence was given. The 
cEieria used by CHOICE in ranking the various channel subsets were described 
Computer time and core allocation were discussed. Finally, some of the ways 
the program has been used at ERL were given. The appendix gives an example 
of a CHOICE computer run with real data. It also provides a description of 
the software including logic flow and Fortran listings. 
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Examples of Output 

Task: Select the "best" 4 of 12 channels for 9 classes, 
(data is from flight line C-1) 

Computation time: 35 seconds 
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Criterion 1 


Maximizing the arithmetic mean of the 
pairwise divergences, (top two sets). 
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ei8f>LAV OF CHANTCLS RANKED ACCCROINC TO AVERACC PAtR~WISE OIVERCENCE. 


CMA»*CLS= t 9 11 12 

AVCRACE PAIR-WISE 01VERGCNCC= 4AS.O 

MINXKUM PAIR-WtSC OIVERC€NCE= 25, « 

RATIO CF THIS OW»CL SET WITH CHAWCL SET YIE1J3IN0 MAxIHtJH PERCENT SEPARATION^ .33 
AVERACC DIVERGENCE BY CLASS... 

AVERACC IKTERCLASS DIVERGENCE FCfi CLASS S = 173.702 

AVERAGE IKTERCLASS DIVERGENCE FOR CLASS C = PCS.OTS 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS O = 180.312 

AVERAGE IKTERCLASS CIVERCCNCE FOR CLASS W = 303.635 

AVERACC IKTERCLASS DIVERGENCE FOR CLASS R = 654.247 

AVERAGE IKTERCLASS CIVEKCENCE FOR CLASS A = 696.084 

AVERAGE IKTERCLASS DIVERGENCE FiOR CLASS Y = 206.235 

AVERAGE IKTERCLASS DIVERGENCE FCR CLASS X = 1165. 065 

AVERACC IKTERCLASS DIVERGENCE FCR CLASS E = 339.305 


PAIR-WISE DIVERGENCE 


CS3 

26.1 06= 

76.4 OC= 

76.6 VS= 

202.4 WC= 

324.9 WO 

110.9 RS= 

209.5 RC= 

100.3 RO 

90.3 RVF 

674.6 

AS= 

186.7 AC= 

82.2 ACc 

92.2 AW= 

733.0 AR= 

25.8 YS= 

117.5 YC= 

204.3 TO: 

37-9 YW= 

52.6 VR= 

423.3 

YA= 

439.3 XS= 

346.8 XC= 

921.4 

804.1 Xthp: 

299.6 XR= 

2962.8 XA= 

3281.7 Xy= 

325.5 ES= 

220.2 EC= 

409.2 

EX^ 

152.0 EU= 

31.5 ER= 

747.5 EA= 

725.9 EY= 

49.6 Cx= 

378.6 






CHAMCLS= 6 9 It 12 

A^RACE PAIR-WlSE DIVERGENCE= 429.7 

MlNtHUH PAIR-WISE 01VERCCKCE= 24.6 

RATIO OF THIS CHAFTCL SET WITH CHAFfJEL SET YIELDING MAXIMUM PERCENT SEPARATION^ .38 
AVERAGE DIVERGENCE BY CLASS... 

AVERAGE IKTERCLASS DIVERGENCE FOR CLASS S = 157.740 

AVCRACE IKTERCLASS DIVERGENCE FOR CLASS C = 259.006 

AVERAGE IKTERCLASS DIVERGENCE FOR CLASS O = 171.964 

AVERAGE IKTERCLASS DIVERGENCE FCR CLASS W = 330.007 

AVERACC IKTERCLASS OIVEHCCNCE FCfi CLASS R = 647.225 

AVERAGE IKTERCLASS DIVERGENCE FCfi CLASS A = 689.260 

AVERAGE IKTERCLASS OIVEfiGCNCE FCfi CLASS Y - 199,601 

AVtKA&C WiTERCLASS DIVERGENCE FCR CLASS X = 1087.277 

AVERAGE Interclass divergence fcr class e = 32s. 354 

PAIR-vasE DIVERGENCE 



27.4 06= 

64.9 OC= 

79.5 WS= 

188.8 WC= 

391.3 WCF 

121.4 

RS= 

225.2 RC= 

82.4 RCB 

96.0 fivp: 

700.0 

AS= 

219.7 AC= 

75.0 AO= 

96.4 AW= 

752.3 AR= 

24.6 YS= 

50.1 

YC= 

133.0 YCe 

67.5 YVf: 

209.0 YR= 

439.8 

TAs 

414.2 X5= 

340.3 XC= 

929.8 XOi 

701.0 XU= 

190.2 XR= 

2831.5 

XA= 

3190.8 XY= 

233.4 ES= 

145.6 EC= 

397.7 

COB 

192.8 EUs 

90.9 ER= 

778.3 EA= 

741.0 ET= 

55.4 ex= 

281.1 
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Criterion 2 


Maximizing the minimum pairwise divergence, 
(top two sets) 





OtaPUAV OF CHANNEL COMBXNAHONS RANKED ACCORDING TO MINIMIS PAIR-WISE DIVERGENCE 


04A»*CLS=: t « to tt 

AVERAGE PAIR-WISE DIVERGENCE^ 304.4 
MIN1»«A4 PAIR-WISE DIVERGENCE^ 

RATIO CF THIS CHAJA^EL SET WITH CMAI»CL SET YIE1.DIN& MAXIMUM PERCENT SEFARATl04= .22 


AVERAGE DIVERGENCE Bt CLASS... 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS S = 
AVERAGE IWTERCLASS DIVERGENCE FCR CLASS C = 
average IMTERCLASS divergence FCR CLASS O = 
AVERAGE INTERCLASS DIVERGENCE FCR CLASS W = 
AVERAGE INTERCLASS DIVERGENCE FOR CLASS R = 
average INTERCLASS DIVERGENCE FOR CLASS A = 
average IKTERCLASS DIVERGENCE FCR CLASS V = 
AVERAGE INTERCLASS DIVERGENCE FCR CLASS X = 
average INTERCLASS DIVERGENCE FOR CLASS E = 


150.230 
231 .255 
140.843 
314.447 
391.805 
423.808 
190.141 
576.936 
320.109 


pair-wise DIVERGENCE 

C9= 35.6 06= 93.1 OC= 122.3 VS= 


AS= 

171.6 

AC= 

iii.e 

AO= 

110.6 

AW= 

YA= 

310.0 

XSs 

152.0 

XC= 

455.1 

W3<= 


143.9 

Etf= 

86.5 

ER= 

609.2 EA= 


231.8 WC= 

665.9 ARs 

416.6 XW= 

64S.S er= 


367.2 WCX= 

38.2 YS= 

331.6 XR= 

57.7 EX= 


89.1 RS= 

121.1 YC= 
1233.8 XA= 

419.5 


184.2 RC= 
21 5'. 0 YO= 

1333.9 XY= 


139.7 RO= 
€3.3 YW= 

272.9 ES= 


87.8 RW= 

181.5 YR= 
212.4 EC= 


0MMC1.S= 1 6 ID 12 

average PAIR-WtSE DIVERGENCE* 295.2 

NINIHLM PAIR-VISE DIVERGENCE* 34.1 

RATIO CF THIS CMAF*^ SET WITH CHAM®. SET YIELDING MAXlKRN PERCENT SEFARATlCW= .15 


AVGRACE DIVERGENCE BY CLASS... 

AVERAGE IKTERCLASS DIVERGENCE FOR CLASS S = 
average IKTERCLASS DIVERGENCE FCR CLASS C = 
average IKTERCLASS DIVERGENCE FOR CLASS O = 
Average ikterclass divergence for class w •= 
average ikterclass divergence FCR CLASS R = 
average ikterclass divergence for class a = 

AVERAGE IMTERCLASS DIVERGENCE FOR CLASS Y = 
AVERAGE IKTERCLASS DIVERGENCE FCR CLASS X = 
AVERAGE ikterclass DIVERGENCE FCR CLASS E - 


152.510 

221.552 

133.666 

305.257 

395.266 

400.699 

196.005 

546.560 

303.221 


pAlR-WtsC DIVERGENCE 
CSS S5-0 06= 89. F 


AS= 

166.5 

AC= 

125-1 

AO= 

YA* 

330.7 

xs= 

S30.6 

XC= 

ro= 

156.2 

EW= 

101.1 

CR= 



122.3 WS= 
99.2 AW= 

394,9 XC= 

571.2 CA= 


228.1 WC= 

612.5 AR= 

582.0 XW= 

563.7 Cr= 


364.9 W5= 

34.1 YS= 

338.3 XR= 

60.6 EX= 


86.6 RS= 

124.3 YC= 

1229.5 XA= 

491.0 


215,6 RC= 

211.6 YO= 

1232.0 XY= 


157.2 EO= 

69.1 TVS 

263.9 CS= 


82.4 RW= 

181.4 TR= 

208.1 EC= 


541.8 
299. 6 
383.3 


529.6 

342.4 

361.9 
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Criterion 3: Maximizing the average divergence by class 

1. Summary print. For example, channel set 
(1, 8, 11, 12) yields the largest average 
divergence for class Rye (Y)-223. Also, 
if this set is used. Bare Soil (X) has an 
average divergence of 924 from all other 
classes. 

2. Expanded print. The first set is (1, 8, 11, 12); 
the "best" set for detecting Rye. The second 
set is (1, 9, 11, 12); the "best" set for 
detecting bare soil. 


OJ 



MATCIX OF 0XVER6ENCES FOR CMANNEt $ET tlELOINO THE MAXIKUH SEPARATION CF A CLASS FROM ALL OTHER CLASSES 



FEATURES 

S 

c 

O 

u 

R 

A 

Y 

X 

C 

t 


11 

le 

S 

174, 

26S. 

180. 

304. 

694. 

€96. 

206. 

1169. 

339 

t 


11 

12 

c 

174. 

26S. 

160. 

304. 

6S4. 

696. 

206. 

1169. 

339 

t 


11 

12 

o 

174. 

26S. 

160. 

304. 

654. 

696. 


1169. 

339 

1 


« 

11 

w 

163. 

247, 

192. 

346. 

489. 

918. 

204. 

743. 

347 

% 


11 

12 

R 

174. 

268. 

160. 

304. 

694. 

696. 

206. 

1169. 

339 

9 


11 

12 

A 

127. 

228. 

197. 

243. 

649. 

706. 

182. 

1133. 

323 

1 


11 

12 

V 

169. 

239. 

166. 

281. 

399. 

607. 

223. 

924. 

363 

1 


11 

12 

X 

174. 

268. 

180. 

304. 

694. 

696. 

206. 

1169. 

339 

1 


9 

11 

E 

16Q. 

244. 

146. 

287. 

489. 

917. 

203. 

751, 

367, 

DIACONAL 



174. 

268. 

ICO. 

346. 

654. 

7D6. 

223. 

1165. 

367. 
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CHANNELSs t « 11 12 

AVERA6C PAIR-WISE DIVERCCNCC= 3dl.4 

MINIlftPI PAIR-WISE OIVERCCNCe= 23.7 

RATIO OP THIS CHAHNCL SET WITH CHANNEL SET YlELDlNO MAXIMUM PERCENT SCPARATION= .24 
AVERACC OIVCRCENCc'BV CLASS... 

AVERACE imERCLASS OlVERCENCE FC« CLASS S = 164.967 

AVERAGE INTERCLASS DIVERGENCE F<::« CLASS C = 233.160 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS O s 165.692 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS W = 281.203 

AVERAGE INTEHCLASS DIVERGENCE FCR CLASS R = 558.932 

AVERAGE XNTERCLASS DIVERGENCE FCR CLASS A = 606.938 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS Y = 222.906 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS X s 923.967 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS E = 362.896 


PAIR-WISE divergence 


cs= 

23.7 OS= 

87.7 OC= 

65.7 W5= 

172.1 WC= 

270.6 VO= 

96.3 RS= 

207.2 RC= 

100.4 RO= 

76.x RW= 

612.7 

AS= 

193.7 AC= 

65.3 AO= 

83.6 AkP 

698.2 AR= 

24.9 YS= 

156.0 YC= 

248.8 YO= 

52.6 YVf= 

50.5 YR= 

422.1 

TA= 

458,6 XS= 

238.3 XC= 

620.5 XO= 

660.0 XkP 

295-3 XK= 

2301.1 XA= 

2521.6 XY= 

344.0 es= 

241.0 CC= 

446.2 

»>= 

183.7 EWP 

54,0 ER= 

727.1 EA= 

769.6 EY= 

50.7 EX= 

410.9 






CHAMCT..S:: 1 9 11 12 

/.VERACE PAIR-WISE 01 '/ERCE^fCE- 443.0 

MINIMUM PAIR-WISE D1VERCENCE= 25.8 

RATIO OF THIS CNAMEL SET WITH CHAI'MEL SET YIELDING MAXIMUM PERCENT S£PARATla't= .33 
AVERAGE DIVERGENCE BY CLASS- . . 

AVERACC INTERCLASS DIVERGENCE FCR CLASS S = 173.702 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS C = 266.075 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS O = 180.312 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS W = 303.633 

AVERAGE INTERCLASS CIVEKCENCE FOR CLASS R = 654.247 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS A = 696.034 

AVERACC INTERCLASS DIVERGENCE FCR CLASS V = 206.235 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS X = 1165.065 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS E = 339.305 


PAIR-WISE DIVERGENCE 


CS= 

26.1 06= 

78.4 cc= 

76.6 WS= 

202.4 WC= 

324.5 WCc 

110.9 

RS= 

209.5 RC= 

100.3 KC*= 

90,3 RW= 

674.6 

AS= 

108.7 AC= 

62.2 AO= 

92.2 AVA= 

733.0 AR= 

25.8 VS= 

117-5 

YC= 

204.3 YO= 

37.9 YVfc 

52.6 YR= 

423.3 

YAs 

439.3 XS= 

346. B XC= 

921.4 XCt: 

8l>4.1 XW= 

299.6 Xk= 

2962.8 

xA= 

3281./ XT= 

025.3 Es= 

220.2 CC= 

4U9.2 

was 

ise.o evp: 

31.3 ER= 

747,5 EA= 

725.9 EY= 

49.6 EX= 

376.6 
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C.’iterion 4 


Maximizing the geometric mean of the pairwise divergences. 
This criterion is also referred to as "per cent separation". 





OISPLAV OF CHANNEL C0»«INAT10NS RANKED ACC0R01N& TO MAXIMUM PERCENT SEPARATION. 


CMAr*CLS= 1 e 9 12 

AVERACE PAIR-WISE DIVERGENCES 393.5 

MINIMUM PAIRHiilSC DIVERGENCES 27.9 

RATIO OF THIS CMAMCL SET WITH CHAf*0. SET YIELDING MAXlWX PERCENT SEPARATIONS l.DO 
AVERACE DIVERGENCE BY CLASS... 

AVERAGE IKTERCLASS DIVERGENCE FCR CLASS S = 165. 56S 

AVERACE INTERCLASS DIVERGENCE FOR CLASS C = 241 .099 

AVERAGE INTEHCLASS DIVERGENCE FOR CLASS O = 146.228 

AVERAGE IMTERCLASS DIVERGENCE FOR CLASS W = 343.142 

AVERAGE IMTERCLASS DIVERGENCE FOR CLASS R = 506.152 

AVERACE INTERCLASS DIVERGENCE FOR CLASS A = 508.696 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS Y = 208.093 

AVERAGE IMTERCLASS DIVERGENCE FCR CLASS X = 760.626 

AVO?ACe IMTERCLASS DIVERGENCE FCR CLASS E = 317.945 


PAIR-WISE DIVERGENCE 


cs= 

27-9 OS= 

74.6 OC= 

86.4 WS= 

246.1 WC= 

404.3 VO= 

111.7 RS= 

214.0 RC= 

109.2 RO= 

88.2 RW= 

646.8 

AS= 

167.4 AC= 

76.4 AC^ 

103.3 AW= 

742.2 AR= 

37.8 VS= 

120.2 YC= 

196.6 YO= 

71.7 YVP= 

204.2 YR= 

396,7 

VA= 

3S6.T XS= 

244.3 XC= 

650.9 XCC 

502.3 XWP 

290.4 XR= 

1845.9 XA= 

1956.9 XY= 

262.2 ES= 

210.1 EC= 

377.1 

03= 

l«/.6 EW= 

99.4 ER= 

710.6 EA= 

608.7 EY= 

56.4 EX= 

333.7 






CHAI«CLS= 1 6 9 11 

AVERAGE PAIR-WISE DIVERGENCES 336.5 

MINIMUM PAIR-WISE DIVERGENCES 27,6 

RATIO OF THIS CMA^^<EL SET WITH CHAhtJEL SET YIELDING MAXIMLM PERCENT SEPARATIONS .89 
AVERAGE DIVERGENCE BY CLASS... 

AVERACE IMTERCLASS DIVERGENCE FCR CLASS S = 162.926 

AVERACE IMTERCLASS DIVERGENCE FCR CLASS C = 246,966 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS O = 151.837 

AVERACE INTERCLASS DIVERGENCE FO? CLASS W = 346.474 

AVERAGE IMTERCLASS DIVERGENCE FOR CLASS R = 489.120 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS A = 517,510 

AVERAGE INTERCLASS DIVERGENCE FOR CLASS Y = 203.546 

AVERACE IMTERCLASS DIVERGENCE FOR CLASS X = 742,910 

AVERACE IMTERCLASS DIVERGENCE FOR CLASS E = 347.346 


PAIR-WISE DIVERGENCE 


cs= 

27.6 06= 

81,2 OC= 

91.2 US= 

249.8 WC= 

420.7 WCb 

110,3 RS= 

187.0 RC= 

90.6 RO= 

96.2 RW= 

655.7 

AS2 

171.7 AC= 

69.0 AO= 

112,3 AW= 

772.5 AR= 

38.1 YS= 

117,5 YC= 

203.0 Y(3: 

64.8 TW= 

196.7 YR= 

364.3 

TA= 

397.0 XS= 

249.1 XC= 

654.2 XO= 

902.4 XW= 

285.7 XR= 

1721,5 XA= 

1897.6 XY= 

272.0 ES= 

223.4 EC= 

419.4 

co= 

196.2 EW= 

60.4 CR= 

759.6 EA= 

721.9 EY= 

53.1 EX= 

364.6 
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Criterion 5: Maximizing the pairwise divergence for each class pair. 


ZO 





0I8PLAV OF CHANNEL COHBtNATlONS TICLOINC MAXlMLW PAIR-WISE OIVERftCMCE 


CLASS MAX A VC 

PAIR DIJ OtJ OWMELS 


cs 



s«.s. 

332.2 

e 

9 

lO 

It 

•s.t 

304.4 

1 

« 

10 

11 

ttt.S 

304.4 

1 

6 

lO 

It 

Zes.s 

242.2 

1 

6 

9 

to 

AAl.Q 

332.2 

6 

9 

to 

11 

12t .4 

429.7 

6 

9 

11 

12 

>2«.e 

282.9 

1 

6 

11 

12 

isr«.2 

295.2 

1 

6 

lO 

12 

«0«,7 

332.2 

6 

9 

lO 

U 

TOO.O 

429.7 

6 

9 

11 

12 

st».r 

429,7 

6 

9 

11 

12 

125.1 

295.2 

1 

6 

lO 

12 

115.2 

228.2 

2 

6 

9 

to 

778,4 

337.3 

2 

6 

9 

It 

38.2 

304.4 

1 

6 

lO 

11 

173.2 

160.5 

1 

2 

7 

6 

273.8 

177.2 

1 

2 

8 

lO 

71.7 

355.5 

1 

6 

9 

12 

233.8 

202.5 

6 

7 

9 

to 

439.8 

429.7 

6 

9 

11 

12 

438.6 

391.4 

1 

8 

It 

12 

361.3 

416.6 

8 

9 

11 

12 

983.3 

403.7 

9 

10 

11 

12 

604.1 

443.0 

1 

9 

11 

12 

338.3 

295.2 

1 

6 

lO 

12 

3008.5 

416.6 

8 

9 

11 

12 

3378.3 

416.6 

8 

9 

11 

12 

544.0 

391 .4 

1 

6 

11 

12 

265.2 

307.4 

1 

2 

8 

11 

469.5 

351.4 

1 

8 

9 

11 

183.7 

391.4 

1 

6 

11 

12 

105.1 

264.4 

6 

r 

lO 

12 

818.2 

332.9 

2 

7 

9 

11 

826.5 

333.3' 

2 

8 

9 

11 

62.4 

282.1 

2 

6 

10 

12 

453.2 

319.1 

1 

0 

ID 

11 
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Show Requests 



Olt^tAY OF CHANNEL COMlINATtCNS RCauesTCO BY *SHCW«* CARDS 


CMAM4CLS9 t 9 tt 12 

AVCRACC PAlR-WlSe D1VERCCNCE= 443,0 

Minima pair-wise oiveRCENCE= as.e 

RATIO OF THIS CHAKHEL SET WITH CHAM«L SET YIELDING MAXIMUM PERCENT SEPARATION^ .53 


AVERAGE DIVERGENCE BY CLASS... 
AVERAGE IKTERCLASS CIVEKCCNCE FCfi 
AVERAGE IWTERCLASS DIVERGENCE FOR 
AVERAGE INTERCLASS DIVERGENCE FOR 
AVERAGE INTERCLASS DIVERGENCE FOR 
AVERAGE INTERCLASS CIVERCENCE FCR 
AVERAGE INIERCLASS DIVERGENCE FCR 
AVERAGE INTERCLASS DIVERGENCE FCR 
AVERAGE INTERCLASS DIVERGENCE FCR 
AVERAGE INIERCLASS DIVERGENCE FCR 


CLASS S = 
CLASS C = 
CLASS O = 
CLASS W = 
CLASS R = 
CLASS A i 
CLASS y - 
CLASS X = 
CLASS E - 


175.702 

263,075 

160.512 

505.635 

654.247 

696.094 

206.235 

1165.065 

559.505 


PATR-WISE DIVERGENCE 


CSs 

26.1 06= 

7B.4 cx:= 

76.6 WS= 

202.4 WC= 

324.5 VQ= 

110.9 RS= 

209.5 RC= 

100.3 

90.3 RWP: 

674.6 

AS= 

183.7 AC= 

82.2 AO= 

92.2 AVP= 

733.0 AR= 

25.3 YS= 

117.5 YC= 

204.5 YO= 

37.9 YW= 

52.6 YR= 

423.3 

VA= 

459.3 XS= 

546.8 XC= 

921.4 »CC 

604.1 XW= 

299.6 XR= 

2962.6 XA= 

3281.7 XY= 

325.5 ES= 

220.2 ec= 

409.2 

to= 

1S2.0 EW= 

51.5 ER= 

747.3 EA= 

725.9 EY= 

49.6 CX= 

373.6 






CHANNEL&S 6 9 lO 11 

AVERAGE PAIR-WISE 01VERCENCE= 532.2 


MINIMIA4 PAIR-WISE DIVERGENCES 26.2 

RATIO CF THIS CMAM'CL SET WITH OWNEL SET YIELDING .“aXIMUM PERCENT SEPARATIONS .04 


AVERAGE DIVERGENCE BY CLASS... 

AVERAGE INTERCLASS DIVERGENCE FCR CLASS 5 = 
average INIERCLASS DIVERGENCE FCR CLASS C = 
AVERAGE INTERCLASS DIVERGENCE FOR CLASS O = 
AVERAGE INTERCLASS DIVERGET4CE FOR CLASS W s 
AVERAGE INTERCLASS DIVERGENCE FOR CLASS R = 
AVERAGE INTEHCLASS DIVERGENCE FCR CLASS A = 
average INTERCUASS CIVERCENCE FOR CLASS Y = 
AVERAGE INTERCLASS CIVESGENCE FOP CLASS X = 
AVERAGE INTERCLASS DIVERGENCE FOR CLASS E = 


136.941 
254.047 
133.203 
332.184 
471 .094 
503.224 
173.453 
653. 091 
332.627 


PAIR-WlsE DIVERGENCE 


CSs 

36. S 06= 

61.5 OC= 

105.8 WS= 

208.3 

1 WC= 

AS= 

169.3 AC= 

102.6 AOS 

109.9 AWs 

753.5 

AR= 

YA= 

339.7 XS= 

233.1 XC= 

675.8 XO= 

378.7 

XW= 

EOS 

131.9 CWs 

73-2 CR= 

770.6 EA= 

739.6 

EY= 



441.0 WO= 

112.4 RS= 

167.0 RC= 

104.3 RO= 

106.7 RW= 

26.2 YS= 

53,7 YC= 

164.4 VOS 

58.9 YW= 

2D6.6 YR= 

172.9 XR= 

1550.2 XA= 

1764.3 XY= 

149.6 es= 

165.8 CC= 

55.6 CX= 

299-7 





664.0 

559.5 

401.2 
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Software Description 


3o 



Overvi ew 


Routi ne 

Function 


CHOICE 

Main program. Coordinates subroutine 

calls . 

INPUT 

Inputs and lists data cards. 


DIVERG 

Computes divergences, does exhaustive 
and saves top-ranked sets. 

search. 

DISPLA 

Displays top channel sets. 


DIV 

Recomputes pairwise divergences for those 
top-ranked sets to be displayed. 



LOGIC FLOW 



32 










Uni vac 1108 Exec 8 Stream 


@RUN,P 

RUN I D , ACCOUNTNO .PROJECT I D , 3 , 1 00 

@ASG,T 

CHOICE,!, 11848 

. MAGNETIC TAPE CONTAINING PROGRAM FILE 

OCOPIN 

CHOICE. 

. TRANSFER PROGRAM FILE TO DRUM 

OFREE 

CHOICE. 

. RELEASE TAPE 

eXQT 

TL0003 

. EXECUTE CHOICE'S ABSOLUTE ELEMENT 

data 

cards 

See: Example of Input Data 


OFIN 


. RUN COMPLETED 



Fortran Listing 


3^ 



CHOICE 





rOA.8 CHOICE, CHOICE 
Eon »4l.-0«/12-0e’42 <2,0> 


MAIN PROGRAM 

«ORACE USED- COCEO> 000026 DATA (0| OQQOll BLANK CC«MCN{2) 000000 
COMMON BLOCKS* 


0003 

IMITS 

000002 

0004 

MAX 

000005 

DDOS 


025416 

0006 

INV 

005514 

0007 

YOU 

015735 


CXIEKNAL REFERENCES (BUXK, NAfC) 


ooto 

IM’UT 

OOll 

DIVERG 

0012 

DISPLA 

0013 

NINTRS 

0014 

NL-CUS 

0013 

NIOIS 

0016 

NtCP$ 

0017 

NSTOPi 


(BUDCK, TTFE, RELATIVE UXATiCN, NAKE) 


STORAGE ASSIGNKENT 

OOOO 000002 tOOT 
0004 OOQOOl MAXCLS 

0004 OOOOOO MXK3MT 

0005 024245 NSHCW 

0007 015733 2AVG 

0007 004000 ZMADIJ 

0007 004050 ZM^C^t J 


OOOl 

OOdOlO 

153G 

0004 

000002 

}«eEST 

0005 

024202 

KBEST 

0003 ] 

t OOOOOl 

PRINT 

0007 

015672 

ZCAVG 

0007 

D03700 

2MCDIJ 

0007 

004120 

ZMPDIJ 


0003 I OOOOOO CARD 

0005 024244 K3I J 

CX305 024201 NCL 

0005 024246 SHZW 

0005 001400 ZCOV 

0007 000760 ZM3IJ 

0007 015732 ZFCT 


0005 

024203 

CHAR 

0004 

ocxxm 

K=RNT 

DQ05 

024200 

W 

0006 

0055DO 

VEC 

0007 

OOOOOO 

2DIJ 

13005 

QQQoaO 

ZfCAN 

0006 

DD47O0 

zv 


0005 024436 01 J 

0004 000003 MSHO 

0005 024243 NPRT 

OOOO 1 DQ0C500 VERS 

OOOS OOOOOO ZINV 

0007 015734 2MIN 


00100 

00101 
f»101 

00103 

00104 
CO104 

00105 

00106 
00107 
DOllO 

00111 

00112 
00113 

00113 

00114 


1« c 

it* 

2* C 
2* 

2* 

» C 

e* 

z* 

2 * 

2* 

Z* 

Z* 

2* 

Z* 

2* 


C. XNE5 NASA/ERL 

INCLUDE SPEC, LIST PARAMETERS, D1^C^KICNS.CC^!HCM 

JCNES 

XW^LICIT INTEGER <A-T) ALL INTEGER EXCEPT Z 
parameter MAXCLA=32 MAX- ND. CT CLASSES. CAN 

„i,-ce-«ci- «ir MEETED. 

“""pARAi^TER MXMSMA^24 MAX ND- OF MEASUtEKENTS (CHATTCtS) . CAN INCREASE. 
PARAMETER MBESTA=12 MAX- HO. OF HEASCREHENTS IN 'BEST* SLOSET 
parameter RSUMA= {M3ESTA+1J4WCESTA/2 StW CF 1,2,..., NT, IN ‘BEST SL©SET 

PARAMETER StWA= (MXMSHA+1 > 4MXMSMA/2 SIH OF 1,2 MAX ND. MEASUREMENTS 

COMHON/UNXTS/CARD, PRINT I/O WIT5-CARD READER, PRINTER 
COMMON /MAX/MXMSMT,MAXCLS,>©EST,MSH0,FFRNT MAXlMlMS 
CCHMCN /INP/ ZHEAN,ZCOV,NM,NCL, 

♦NBEST, CHAR, NPRT. MDrj.NSTW.SPOJ.CIJ INFin- 

DIMENSION 2M£AN(MXMSHA,MAXCLA) HOLDS MEANS FOR EACH MEASWieHENT, CLASS 



i6 



OOtl» 

OOllS 

OOtlG 

cxniT 

ootza 

00t2t 

00121 

00122 

00122 

00123 

00124 

00125 
0012ft 

00126 
00127 

00130 

00131 

00132 

00133 

00134 

00135 

00136 
00157 

00140 

00141 
D0141 

00142 

00143 
00143 
00147 
00151 
00157 
00160 
ooieo 
00161 
00161 
00162 
.00163 
00164 


26 

2* 

2* 

e« 

2* 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

26 

36 

46 

36 

66 

76 

66 

96 

106 

116 

126 

136 

146 

156 


DIlCNStQN ZCOVtSUMA.MAXCLA) HCL03 COVARIANCES (LOUER TRIANOCE) 

C FOR EACH CLASS. 

PARAMETER PCTPRT=10 NO, OF COMBO’S RANKED BY PERCENT TO BE PRINTED 
DtHENSiCN VEC(I^ESTA) HOLDS UNIQUE CC6S1NAT10N CF FEATURES (CHANNELS) 
COMMON /INV/ZINV.ZV.VEC INVERSION ARRAYS 

OIHENSION ZINV(RSUHA.MAXCLA) HOLDS INVERSE COVARIANCE MATRICES 
C FOR EACH CLASS 

Ot»CNStON ZV<MBESTA,MAXCLA) VORKINO ARRAY USED VJKEN CCMrUTlNO 
C INVERSE By BCRCERINC METHOD. 

PARAMETER OSUM= (MAXCLA-D6MAXCLA/2 ’MAXCLA* COMBO’S TAKEN TWO AT A TIME 
DIMENSION ZDIJ (C5UM) WORKING BUFFER FOR INTERCLASS DIVERGENCE 
DIMENSION CHAR (MAXCLA) MCLOS CHARACTER REPRESENTATION CF CLASSES 
C06CN /YOU/ ZCIJ.ZMDIJ.ZMCDIJ.ZMADIJ, 

62M)CIJ,ZMFD1J.2CAVC,2PCT,2AVG,ZMIN FOR CIJ 
PARAMETER MSHQUPilO MAXIMUM NO. CF SrtX REQUESTS HONORED 
0I»ENS1CN SHOW ()«ESTA,K5M0W) HOLDS CHAt^fd. 00»®0"S REQUESTED 
PARA>eTER MPRINT=20 MAXIMUM MO. CF PRINT REQUESTS TO BE HONORED 
DIFCNSION 2>01J(DSUM,3) HOLDS MAX. DIJ.AVG. DIJ, PACKED CHAMd. U::RD(FCW) 
OIMENSICN ZMCDl J (MAXCLA, 2) HOLDS MAX. CIJ FOR CLASS |PCW 
OIPCNSICM ZMADIJ (KFRIMT.e) HOLDS MAX, AVO. CIJ.PCWFCfi MO. PRINT REQUEST 

DIMENSION ZhFCtJ IMPRINT, 2) HOLDS MAX. HIN. CIJ.FCWFOR NO. PRINT REQUEST 

PARA«:tER CSUMl=DSl«*l CLASSES TAKEN 2 AT A TIKE ♦ 1 WORD 
DltENSlCN Zrt=CI J (CSLWl .FCTFRT) DU’S RANKED BY FERCEKT.FCW 
OIKENSICN ZCAVG (MAXCLA) WORK ARRAY USED TO FIMJ LARGEST CLASS AVC. DIJ 

OKCNSICN DU(CSUM) HOLDS CHARACTER OOEO’S FOR 01 J PAIRS 

END 

DNCNSICN vers (2) LAST T1)C ANY CHANGE TO PROGRAM 
DATA VERS/’AFR 4’,’, 1973’/ 
data card/ 5/ CARO READER UNIT 

DATA FRINT/6/ PRINTER UNIT 
VRITE (PRINT, lOO) VERS PRINT VERSION MO. 
lOO FORMAT (*0’,2A6,/, ’O’) 

irvPUT GET f'EANS, COVARIANCES, MO. CF fEASUKEMEI4TS, (CHA>#JELS) , 

C MO. OF CLASSES, OIPFNSICN CF ’BEST" SUBSET, DISPLAY UFO. 

CALL OlVERC CCPFim: divergence FCR ALL CCPBINATICNS Af® SAVE 
C THOSE CF INTEREST 

CALL OISPLA OUTPUT SEFARATICNS 

STOP 

EM3 


EMD OF OCHPILATICN' 


MO DIAGNOSTICS. 



INPUT 



rOR.S INPgTitNPUT 

roft •*L-0«/t2'Oe>44 <2,0) 


MJeROUTINE INPUT ENTRY POINT 000920 


STORACE USED' CCCEtt) 000994 OATACO) OOOZS4 BLANK COM3NC2) OOOQOQ 
OOMHON BLOCKS' 


0009 

IMtTS 

000002 

0004 

MAX 

000009 

0009 


025416 

0006 

INV 

005514 

OQOT 

YOU 

019739 


EXTERNAL REFERENCES <BLOCK, NAME) 


OOlO 

. liRCUS 

OOll 

NICSS 

0012 

IMiUS 

0019 

7CRR6* 

0014 

NIOIS 

0019 

NERR3S 


^liCRACC ASSIOi'WENT <BLOCK, TYFEf RELATIVE LOCATION, NA>C) 


oooo txiocnr lor 

0000 OOOt74 126F 

0001 000109 175C 

OOQl OOOl 97 224^ 

OOQl 000900 909C 

OOOl 00049S 9640 

OOOO ooQori eoF 

oooo 000149 99F 

OOOO 1 000006 I 

0004 1 000002 ^GEST 

0004 t DOOOOO MXKSMT 

0005 I 024245 NSHCW 
OOOO I DOOOOO SUM 
0007 OOOOOO ZOIJ 

0005 R OOOOOO z>e:an 

0006 004700 2V 


OOOQ 


000146 

IQQF 

OOOO 


000200 

135F 

OOOO 


000013 

20F 

OOOl 


000167 

233G 

OOOl 


000353 

331C 

OOOl 


000462 

402G 

OE301 


000240 

SOL 

OOOl 


000310 

SSL 

ooca 


000233 

INJP9 

0005 

t 

024244 

7CIJ 

0005 

X 

024202 

^BEST 

0003 

1 

OOOOOl 

PRINT 

0006 


005500 

VEC 

0006 


ooccoo 

ZINV 

0007 


015734 

ZMIN 


OOOO 


000195 

IlOF 

oooo 


000202 

136F 

OOOl 


0001 lO 

201G 

OOOl 


000176 

2370 

OOOl 


000360 

3350 

0001 


000467 

4060 

oooo 


000075 

SIF 

0003 


OOOOOO 

CARD 

OOOO 


000002 

J 

oooo 


000005 


0005 


024201 

NCL 

oooo 


000004 

PTS 

0007 


015733 

ZAVG 

0007 


004000 

ZMACIJ 

□007 


004050 

ZK>^IJ 


OOOO 000166 119F 

0000 000210 150F 

0001 000134 2116 

OOOl 000065 30L 
OOOl 000412 3500 

0000 000064 50F 

0001 000257 92L 

0005 t 024203 CHAR 
OOOO I 000003 K 

0004 I 000004 MlrRNT 

0005 I 024200 NM 

OOOO I OOOOOl RSUH 
□007 015672 7CAV6 

0007 003700 ZMCDIJ 

0007 004120 ZKPOIJ 


□000 


000172 

125F 

oooo 


000221 

160F 

0001 


□00141 

21 5G 

OOQl 


000273 

3010 

OOOl 


000417 

394C 

oooo 


000066 

70F 

oooo 


000120 

93F 

0009 


024436 

DIJ 

0004 


OOOOOl 

MAXCLS 

0004 


000003 

MSHO 

0005 


024243 

NFRT 

0009 


□24246 

SHOW 

0005 

R 

001400 

zcov 

0007 


000760 

ZMDIJ 

0007 


015732 

H=CT 


OOlOt 

1* 


SUBROUTtt^ I^PUT 

ocnoi 

K* 

C 


CX1103 

3* 


INCLUDE SPEC, LIST 

cx>103 

S* 

C 


pot 04 

3« 


IMPLICIT INTEGER CA-Y) 

poibs 

9* 


PARAMETER MAXCLA=32 

pOSDS 

34 

c 

increase if needed. 


PARAMETERS , 0 1 ^ENS IONS , C0MM3N 

ALL INTEGER EXCEPT Z 
MAX. NO. OF CLASSES. CAM 


JONES 

JCNES 



00106 

3* 


DQlOr 

36 


OOltO 

36 


OOltl 

36 


00112 

36 


00113 

36 


00114 

36 


T»il4 

36 

. . 

00113 

36 


00116 

36 


00116 

36 

C 

00117 

36 


00120 

36 


00121 

36 


00122 

36 


00122 

36 

C 

00123 

36 


00123 

36 

c 

00124 

36 


00125 

36 


00126 

36 


00127 

36 


00127 

36 


00130 

36 


00131 

36 


00132 

36 


00133 

36 


00134 

36 


00133 

36 


OD136 

36 


00137 

36 


00140 

36 


00141 

36 


00142 

36 


00142 

36 

EM? 

00143 

46 


OOlSO 

56 

lO 

00151 

66 


00152 

76 


00153 

86 


00154 

96 


00135 

106 


00156 

116 


00157 

126 


00160 

136 


00162 

146 


00172 

156 

20 

c«ir 2 

166 


□0172 

176 


D01T2 

186 


00173 

196 


00174 

206 

30 

001 77 

216 

40 

00206 

226 

SO 

00207 

236 


00210 

246 


00213 

256 

60 


PARAttTCR MXMSMA=24 HAX MO. OF «TASUR£MENTS (CHANNELS) . CAN INCREASE. 
PARAMETER )«ESTA=12 MAX, NO. CF ttASlREHENTS IN ‘BEST’ SlBSET 
parameter RSUMA=tMSeSTA*n4HBESTA/2 SLW CF 1,2,..., NO. IN *OEST- SUBSET 

PARAMETER SLIKA= (MXMSMA»l)*HxHSHA/2 SIW CF MAX NO. MEASUREMENTS 

COMMC«/UNITS/CARD, PRINT I/O UNITS-CARC READER, PRINTER 
COMMON /MAX/MXMSHT,KAXCLS,MBEST,HSH0,W=RNT MAXlMl^S 
COMMON /INF/ 2KEAN,ZCOV,NH,NCL, 

*»®EST,CHAR,M=RT.K>lJtNSH3W,SHCV,ClJ 1M=UT 

OlKENSlON ZKRANIMXMSMA.MAXCLA) HCCCS )EANS FOR EACH MEASUREMENT, CLASS 
DIMENSION ZCOV(S(JMA,MAXCLA) HCtDS COVARIANCES (LCWER TRIAN&LE) 

FCR EACH CLASS. 

PARAICTER PCTFRT=10 NO. CF CC»eO*S RAMUED BY PERCENT TO BE PRINTED 
DIMENSION VEC(PBESTA) HOLDS LNIODE COBINATION CF FEATURES (CHAf-ffCLS) 
COMMON /INV/ZINV,2V,VEC INVERSION ARRAYS 

OIPENSION Z1NV(RSUMA,MAXCLA> HOLDS INVERSE COVARIANCE MATRICES 

FOR EACH CLASS 

OHCNSICN ZV()©ESTA,MAXCLA) WORKING ARRAY USED VNEN CCMPLJTING 

INVERSE BY BORDERING PETHOD. 

PARAFETER CSCK^tMAXCLA-l)4NAXCLA/2 ’KAXCLA’ COEO'S TAKEN TVO AT A TI»€ 
OUCNSION ZDIJ(DStJM) VCfiKING BUFFER FCR INTERCLASS DIVERGENCE 
OIPENSICN CHAR(MAXCLA) HOLDS CHARACTER REFRESENTATICN CF CLASSES 
CCK)CN /YCIJ/ ZCI J,Z)^I J iZMCCIJ.ZKACIJ, 

♦ZMFClJ,ZHf=CIJ.ZCAVC,2FCT,ZAVG,ZMIN FCR DIJ 
PARAMETER MSHCW=10 HAXlM>« NO. CF SHOW REQUESTS HCNRED 
DIFENSlCN SHOW (NBESTA.MSHCW) HXDS CHAhML CCPEO"S REQUESTED 
PARA^ETER M--RINT=2a MAXIMUM NO. CF PRINT REQUESTS TO BE HCNOUED 
OlfENSlCN ZM3U<CSUM,3) HCtDS MAX. DIJ.AVC. DU, PACKED CHA)NEL WDRD(PCW) 
DIPENSION 2HCCIJ(MAXCLA,2> HOLDS MAX. DU FCR CLASS.FCW 

OIPENSICN 2MAD1J(W=RINT.2> HXCS MAX. AVG. D1J,FCW FCR NO. PRINT REQUEST 

OUENSICN 2KFCUtPFRIMT,e> HOLDS MAX. MIN. DIJ.FCW FOR NO. PRINT REQUEST 

PAHAPETER CSUMUOSUJM+1 CLASSES TAKEN 2 AT A T1)E ♦ 1 VERB 
DIPENSICN ZMrDU(CSUMl ,FCTFRT) DU'S RAMUED BY PERCENT. F<V 
OIPENSICN ZCAVG(KAXCLA) WORK ARRAY USED TO FirC LARGEST CL^USS AVG., OU 
ODCNSICN BU(DSUM) HOLDS CHARACTER CCPCO'S FCR DU PAIRS 

read (CARD, 10> PW.NCL.FBEST NO. MEASIEENENTS.MO. CLASSES . SUBitr RAMU 
F0UMAT(SX,I2,SX»12,SX,IZ) CHAKPELS=XX CLASSES=XX best-xx 
KAXCLS=MAXC»_A OCT MAX. CLASSES FRCM PARAPETER STATEMENT 
MXMSMT=MXMS«A DITTO FCR MAX. MO. )EASIREFCNTS 
SUM=SUMA DITTO FOR SUM CF 1,...,MAX. NO. »EASlREMENTS 
KBEST=MBESTA DITTO FCR NO. CF FEANS IN "BEST' 

MSHCFMSHCW DITTO FOR MAX. SHCW REQUESTS 
|ifPl^tslT=MPRINT DITTO FCR HAX. PRINT REQUESTS 

RSUM=RSUMA DITTO FCR SUM CF 1,2,.... no. in -best* .. 

IF(MXMSMT.C€.NM.ANB.MAXCLS.CE.NCL.ArC.FSEST.CE.FBESnGO TO 5Q DIi(l S CK 

WRITE (PRINT, 20) MXMSMT.NM.KAXCLS.NCU.KBEST.NBEST NO, ERROR 
FCRMATC MAXIMUMS EXCEEDED. CHECK CARO CCLUffm^OF I^UT.',/, 

*’ NO. FEASU?;EHErj73 AL!-CWED= " , 1 3 , ' ATiCF-'-.-i' 

NO. CLASSES ALLCVEB=-,U.* ATTEMPTED INPUT=' ,013 ,/, 

♦- max. rank subset ALL0V£C=*.13.- ATTEMPTED INFUT='.013) 

RCTURN 0 ERROR EXIT 

DO 40 J=1,NCL INCEX THRU CLASSES, READ MEANS 

READ (CARO, 30) <2KEAN(K,J) ,K=1 ,NM) READ MEANS 

NO. POINTS IN LCVCR TRIANGLE COV. MATRIX. SUM 1,2,\.,NM 
OO 60 J=1.NCL INDEX THRU CLASSES- READ COVARIANCES. 

READ (CARO. 30) (ZCOV £K. J) ,K=1 ,PTS) READ LCVCR TRIANGLE COV. 
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00215 

26* 

C RE' 

00222 

27* 


00230 

26* 

70 

00251 

29* 


00232 

50* 


00255 

31* 


00230 

52* 


00241 

53* 


00242 

54* 


00245 

35* 

55 

0024« 

56* 


ocess 

57* 

80 

00255 

58* 

C 

00255 

59* 

c 

00255 

40* 

C CH 

00254 

41* 


00256 

42* 


00263 

43* 

91 

00263 

44* 


00264 

45* 


00265 

46* 

90 

00267 

47* 


00274 

48* 

93 

00274 

49* 


00275 

50* 


00276 

51* 

92 

00300 

52* 


00305 

55* 


00311 

54* 

95 

00312 

95* 

94 

CC314 

SC* 

90 

00316 

57* 

lOO 

00316 

58* 

c « 

00317 

59* 


00324 

eo* 

no 

00325 

61* 


00327 

62* 

115 

00330 

63* 


00335 

64# 

120 

00342 

65* 

125 

00543 

66* 


00344 

67* 


00546 

68* 

126 

00547 

69* 


00352 

70* 

150 

00361 

71* 

135 

00362 

72* 


00370 

73* 

136 

00371 

74* 


OD376 

75* 

150 

00377 

76* 


00401 

77* 


00404 

ra* 

155 

00413 

79* 

160 

00414 

so* 


00415 

61# 



READ tCAR0»70) (CHAR (J) , J=1 |NCL) CLASS£S=C ,0» • . . 

FORMAT (rx.3€(lX|Al)) 

INITIALIZE INDEX FOR DU CHARACTER PAIR 
DO 35 1=2, NCL THE 1 OF DU 
R=l-1 THE J IS ALWAYS SMALLER THAN THE I 
DO 55 J=1,R THE J CF DU 

DU(MM)=CHAR(I> GET FIRST CHAR OF CLASS PAIR 

FLDt6,S,0U(MM)l=FLC(O,6,CHAR(J)) GET AM) PACK 2-M> CHAR. IN CLASS FAIR 
PAIR ITCEX 

READ (CARD t SO) I^RTiMDU jNSHCW FRINT=XX,HAX DU=XXXXX,NO. SHCW=XX 
TORMAT(6X,I2.9X,15,I0X,I2) NO. CF COeO'S TO PRINT, MAXIUM INTERCLASS 

DIVERGENCE ALLOuEC, ND. CF SPECIAL SH!M 
CCPeiNATICNS REQUESTED. 

ECK NO. CF PRINT A^C1 SHCVI REQUESTS 
lF(fF«T.LE. W=RKT1 CO TO 90 WD. OF PRINT REQUESTS EXCEEDED MAX. 

(PRINT ,91) f'FRT »NPRT .►FRNT YES, PRINT WARNING 
FCRMATCO MO. CF PRINT RE'.’;.E5TS=* ,15,014,* TOO LARGE, MAXIMUM* 

•• IS=*,13,‘ print F-COUESTS set TO MAX.*) 

M=RT=FrRNT SET PRINT REQUESTS TO MAXIMUM 

IF(NSI-CW.LE.MSHO) CO TO 92 NO. CF SHOW REQUESTS EXCEEDED MAX. 

VRITE (PRINT, 93) NSHOV.NSHJW.MStO YES, PRINT WARNING 
FCRMATCO NO. OF SHCW REQUESTS^* ,15,014, * TOO LARGE, MAXIMUM* 

#* ALLOWED IS=*,I3,“ SHOW REQUESTS SET TO MAX.*) 

NSM0W=MSH0 set show REQUESTS TO MAXIMUM 
IF(NSHCW.EQ.O) CO TO 99 ANY SHCW CARDS 
DO 94 J=1 .NSHOW INDEX TMTU SKW CARDS 

KEAD(CARD,95) (SH0W(1 , J) , 1=1 ,1BEST) OWMELS=Ol ,05,10, .. . 

FORMAT (6X,24 (IX, T2) ) 

CONTINUE 

WTITEtPrCn'.lOC) 

format (IHl,- INTERFRETATICN CF CARD IM’in'*) 

TINT CHANUGLS,CLA55£S.BEST CHATMELS 

VMITE (PR I NT , 1 1'-l > M4 , NCL , fBEST 

FCRMAT(*0 CHA7f'EUS=",I2,* CLASSES=" , 12, * BEST=*,12) 

WRITE (PRINT, 115) 

FCRMATCO IM’OT WEANS*) 

DO 120 J=1 ,NO- INDEX THRU CLASSES 

WRITE (PRINT, 125) (ZMEAN(K,4) ,K=1 ,WM) PRINT MEANS 
FORMAT (5x,5El 5.8) 

PTS=(NM*l)4NM/2 NO. POINTS IN LOWER TRIANGLE COV. MATRIX 
WRITE (PRINT, 126) 

FCRMATCO INPUT COVARIANCE*) 

DO 130 J=1,NCL IfBLX THRU CLASSES 

WRITE (PRINT, 135) (ZCOV <K , J) ,K=1 ,PTS) 

FCRMAT < 5X , 5E1 5 , B) 

WRITE (PRINT, 136) (CHAR (J) , J=1 ,NCL) 
format (*□ INPUT CLASSES ",36(1X,A1)> 

WRITE (PRINT, 1 30) NPRT.MDIJ.NSWCW PRINT=XX MAX. CU=XX NO. SHCW=XX 
FORMAT(“0 PR1NT=*,I2,* MAX CU=*,I5,* NO. SHDW^*,I2) 

IF (NSHOW. EQ.O) RETURN ANY SWCW CARDS 
DO 135 J=l,NSHOwr INDEX THRU SHX CARDS 
WRITE (PRINT, 160) (SMOWd , J) , 1=1 ,NBEST> BEST CHANNELS 
FORMAT (*0 SHCW CHANNELS^* ,24 (IX, 12)) 

RETURN 
CM) 


PRINT LOWER TRIANGLE COV. 
PRINT CLASSES 
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DIVERG 



POA.ft DtVCRC.OIVCRC 
rCR 9«C-06/lZ~0e’44 (1,0] 


SUBROUTIHC OlVEKC CNTRY POINT D0174T 


nORAilC USCO* CODE (11 OQ1T74 DATA(01 000520 OLAM( COMMON (2) 000000 
CO»«ClN BLOCKS’ 


0003 

IWITS 

000002 

0004 

MAX 

000003 

(xns 

INP 

02S416 

cx>os 

INV 

cns314 

0007 

YOU 

015735 


CXTEKNAL REFERENCES (BLOCK, NAME) 


OOlO 

M^CUS 

oaii 

Nices 

0012 

NIOIS 

0013 

tCRRGt 

0014 

TEKR3S 


STOUce ASSICJ^tNT (BLCCK , TYPE, RELATIVE LOCATICN, NATO 


OOOl 


000026 

1466 

OOOl 


(xnii6 

2046 

OOOl 


000211 

2316 

oooo 


000247 

237F 

OOOl 


{XXJ243 

2426 

(xni 


000346 

2716 

OOCXl 


{X»354 

31 IF 

OOOl 


DCX3540 

3506 

OOOl 


00070Z 

4270 

OOOl 


000102 

SOL 

OOOl 


001313 

5756 

(xni 


001404 

6356 

OOOl 


001511 

7066 

odoi 


001665 

7656 

OOOl 


001307 

65L 

OOOl 


001323 

S4L 

(X303 


024436 

OIJ 

OOOO 


000132 

J 

OOOO 


000176 

L 

(XX30 


000201 

m 

0005 


024201 

TiCL 

0003 


OOOOOl 

PRINT 

OOOO 


000162 

VC 

oooo 


000166 

VJl 

oooo 


000140 

XC 

oooo 


000X41 

Xl 

OOOO 


000146 

XK 

OOOO 


000147 



(X»l 


000040 

1546 

OOOl 


000132 

2106 

oooo 


000225 

232F 

oooo 


000256 

238F 

oooo 


OOQ322 

243F 

OOOl 


001146 

300L 

oooo 


000362 

312F 

OOOl 


000554 

3S4C 

OOOl 


000720 

4416 

OOOl 


001123 

5046 

OOOl 


0Q1330 

6046 

OOOl 


0014’- ( 

S426 

OOQl 


001550 

7236 

OOOl 


001666 

7700 

CXXll 


001423 

89L 

OOOl 


000036 

96L 

oooo 

I 

OOOOOO 

ELE 

oooo 

I 

000135 

JJ 

oooo 

I 

000200 

LOW 

0004 


000004 

MFRNT 

GC05 

I 

024200 

fc; 

0005 


024246 

SHOW 

oooo 


000157 

vjc 

oooo 


OOOl 70 

WKL 

oooo 


000136 

XCL 

oooo 


000150 

XU 

oooo 


000155 

XKJ 

oooo 


000156 

XPTS 


OOOl 


000064 

1646 

OOOl 


000160 

2226 

oooo 


000243 

233F 

oooo 


000266 

239F 

oooo 


000327 

244F 

(XXII 


000407 

3046 

(XXX) 


000400 

31 3F 

OOOl 


000600 

3666 

□OOl 


000740 

4456 

OOOl 


□01164 

5326 

□UOl 


001343 

6126 

OOOl 


(XH434 

6556 

OOOl 


001602 

7376 

OOOl 


001701 

7766 

OOOl 


001 465 

SOL 

OOOl 


001707 

99L 

oooo 

1 

000202 

I 

oooo 

I 

000204 

JJJ 

0004 


(XX5001 

HAXCUS 

0004 


000003 

HSHO 

0005 

I 

024243 


oooo 

1 

000203 

UP 

oooo 

I 

D00164 

WI 

oooo 

I 

000167 

WL 

oooo 

I 

000137 

XCT 

oooo 

1 

000144 

XtND 

oooo 

t 

000151 

XKJ 

oooo 

R 

000172 

Z 


OOOl 


000104 

2006 

OOOl 


000334 

203L 

OOOl 


000503 

23QL 

OOOO 


000206 

231F 

OOOl 


CKK1233 

2346 

OOOO 


QOOZ45 

236F 

OOOl 


0G0710 

240L 

oooo 


000303 

241F 

OOQl 


000256 

2456 

OOOl 


000310 

2606 

0001 


000424 

3076 

OOOO 


000331 

310F 

OOOl 


000456 

3176 

OOOl 


000531 

3416 

OOOl 


000070 

4QL 

OOOl 


000660 

4176 

OOOl 


000762 

4536 

OOOl 


000775 

4566 

OOOl 


001243 

5516 

0001 


□01274 

3636 

OOOl 


001356 

6216 

OOOl 


001365 

6256 

OOOl 


001451 

6636 

OOQl 


001474 

7006 

□OOl 


001626 

7506 

0001 


001637 

7546 

OOOl 


001227 

82L 

OOOl 


001231 

83L 

OOOl 


001463 

91L 

OOOl 


001525 

93L 

0003 


OOOOOO 

CARD 

0005 

1 

024203 

CHAR 

OOOO 

I 

000205 

II 

OOOO 


000424 

INJF* 

oooo 

I 

000133 

J2 

oooo 


000134 

J3 

0004 


(moc;r >2 

MBEST 

0005 


024244 

KCIJ 

0004 


OODOOO 

MXMSMT 

0005 


024202 

NBEST 

0003 


024245 

NG}»3W 

OOOO 


000174 

PAIRS 

0006 


005300 

VEC 

oooo 


000160 

WA 

oooo 


000171 

VfU 

oooo 


D0O165 

WJ 

oooo 


000173 

VM 

OOOO 


□00161 

V,(JP 

oooo 


000142 

XD 

OOOO 


000143 

XF 

oooo 


000143 

XJ 

oooo 


000153 

XJl 

oooo 


000154 

XKK 

oooo 


000152 

XL 

0007 

R 

015733 

ZAVC 

0007 

R 

: 015672 

ZCAV6 



OOOS R CU 1400 ZCOV 
DQOT R 004000 ZMACIJ 
0007 R 004090 ZHHDtJ 
oooe R 004700 ZV 


0000 R 000163 ZO 
0007 R 003700 ZHCD14 
□007 R 004120 ZKPOIJ 


0007 R OOOOOO ZCIJ 
0007 R 000760 ZMOtJ 
0007 R 015732 ZFCT 


0006 R OOOOOO ZINV 
0005 R OOOOOO Z^tEAN 
0000 R 0001 75 2PCW 


0000 R 000177 2H 
0007 R 015734 2MtN 
0000 R 000116 ZPRNT 


OOlOl 

1 * 


00101 

24 

c 

DD 103 

34 > 


CX 1103 

44 > 

c 

00103 

44 

c 

00104 

44 


00105 

44 


ootos 

44 

c 

00106 

44 


00107 

44 


00110 

44 


00111 

44 


00112 

44 


00113 

44 


00114 

44 


00114 

44 


00113 

44 


00116 

44 


00116 

44 

c 

00117 

44 


001 eo 

44 


00121 

44 


00122 

44 


00122 

44 . 

c' 

00123 

44 


00123 

44 

c 

00124 

44 


00123 

** 


00126 

44 


00127 

44 


00127 

44 


00130 

44 


00131 

44 


00132 

44 


00133 

44 


00134 

44 


001 3 S 

44 


00136 

44 


00137 

44 


□0140 

44 


00141 

44 


0014 P 

44 


00142 

44 


00143 

94 


00143 

64 

( 

00144 

74 


00144 

64 

< 

00145 

94 


00130 

104 

) 


SUBROUTlhC: OIVERO 

INCLUDE SPEC, LIST PARAMii JERS,D1^0«SICNS,CCMMCN 

INITIALIZE -VEC* FOR UNIRUE INTEOGl CCNSINATICN ALCCRITW 


JONES 


JCMES 


ALL INTEOER except Z 

max. no. of classes, can 


implicit integer (A-Y) 

PARAIEITER MAXCLA=S2 
INCREASE IF ^EEDED. 

parameter MXMSMA=Z4 
PARAIETER M5ESTA=12 
PARAMETER RSUMA= (M2ESTA+1 ) iPBESTA/2 
PARAWETTER SUMA= (MXMSMA+l) 4MXMSMA/2 


MAX r«5. CF lEASUREKENTS (CHAl^ELS) . CAN INCREASE. 
MAX. NO. CF f-EASCnOCNTS IN "BEST* SLBSET 

SUM CF 1 , 2 ,...,N 3 . IN "BEST" SUISET 
SUM CF 1 , 2 ,...»MAX NO. MEASUREMENTS 
CCHM 3 N/LN 1 TS/CARQ, PRINT I/O LNITS-CARD REACERj PRINTER 

OWCN /KAX/MxMSMT,MAXCLS,MBE 5 T.MSN 0 ,MPRNT MAXIMUM 
COPMCN /IhF/ 2 VEAN,ZCOV,NM,NCL» 

««est,char,nfrt,mdij,nshcw,sh:w,dij input 
OIPENSION 2VEAN{MXMSMA,MAXCLA) HCLCS MEANS FCR EACH MEASUREMENT, CLASS 
dimension ZCOYISUMA.MAXCLA) IOXS covariances (LCVER TRIANGLE) 

for each class. 

parameter PCTFRT=10 no. CF CC7©0‘S RAMKED BY PERCENT TO BE PRINTED 
DIMCNSICN VEC (MBESTA) MCLCS UJIQUE CCMBINATICN OF FEATURES <CHAM»ELS) 

0cm34 /INV/ZINV,ZV,VEC INVERSION ARRAYS 

DIMCN31CN ZINV(RSUMA,MAXCLA) HOLDS INVERSE COVARIANCE MVITRICES 

FOR EACH CLASS 

DIMENSION ZVCMBESTA.MAXCLA) VERKING ARRAY USED M«EN 

inverse by bcrcering mcthco. 

PARAMCTER DSUM= IMAXCLA-l ) ««AXCLA/2 "MAXCLA" CCMBO'S TAKEN TWO AT A TltC 
DIMENSION ZOIJ(DSUM) WORKING BUFFER FOR INTERCLASS 
OIMENSICN CHARIMAXCLA) HOLDS CHARACTER REPRESENTATION CF CLASSES 
OCMIMCN /YDIJ/ 2D1J,ZMCIJ.ZMCBIJ,ZMADIJ, 

«MM3IJ,ZMPCIJ,ZCAVG,ZPCT,ZAV0,ZMIN for dij 

parameter MSHOW=10 maximum no. CF SKM REGIESTS HCNCREO 
DIMENSION S«:W <M«ESTA,MSHOW> HCLDS CHAMtf^ CCMfiO-S REOLESTED 
parameter »fi=RINT=20 MAXIMUM NO, CF PRINT REQLESTS TO BE HCNCRED 
DIMENSION ZHDIJ(CSUM,3) HCLCS MAX. DU.AVO. DIJ, PACKED CHANNEL WCRDIPCW) 
DIMENSION ZHCCIJ (MAXCLA. 2) HOLDS MAX. DIJ FOR CLASS, PCW 

D ^ ICN 2MADIJ(MPRINT.2> HOLDS MAX. AVG. CU.FCW FCR NO. PRINT REQUEST 

ZMHDW<^rfiINT.^> HCUCS MAX. H.N. O.A.PCW r<K ™.MI RMUIST 

PARAMETER DSUM1=CSUM«-1 CLASSES TAKEN 2 AT A TIME + ^ VCRD 
OiHEKSiaJ rrt-CI J iCSU-;l .FCTFST) CI J*S ranked by PERC...i7 .rCW 

WMENSION ZCAVG (MAXCLA) WORK ARRAY USED TO FIM3 LARGEST CLASS AVG. DIJ 
OIMENSICN OIJ(DSUM) MOLDS CHARACTER CCMSO'S FOR DIJ FAIRS 


END 


dimension ELE(RSOMA) PRINF BUTER. holds CLEMENTS CF 

COVARIANCE M1ATR1X ON ERROR. 

DIMENSION ZPRNT (M3ESTA) PRINT BIFFER. MCLDS ROW 
OP COVARIANCE MATRIX ON ERROR. 

DO 10 J=1,NBE5T DIMENSION IS "BEST* 

VEC<J)=J SET TO FIRST INIQl^ SET 



00152 

114 


00152 

124 

C 

00153 

134 

96 

00156 

144 


00157 

154 


00160 

164 


D0161 

174 


00163 

1S4 


00166 

194 

30 

OOITO 

204 

40 

00172 

214 

20 

00174 

224 


00175 

234 

50 

00176 

244 


00176 

244 

C 

00176 

244 

c 

00176 

244 

c 

00176 

244 

c 

00176 

244 

c . 

00176 

244 

c 

00176 

244 

c 

00176 

244 

c 


001T6 

001T7 

UOZOZ 


x>cs 


244 

244 

244 


00203 

244 


00206 

244 


00207 

244 


00212 

244 


00213 

244 


00214 

244 


00219 

244 

221 

00219 

244 

C N 

00221 

244 


00224 

244 


00229 

244 


00229 

244 

C 0 

00227 

244 


00230 

244 


00233 

244 


00236 

244 

202 

00240 

244 


DDZ41 

244 


00244 

244 


00247 

244 


00250 

244 


IXJ291 

244 

204 

00293 

244 


00296 

244 

209 

00297 

244 


00262 

244 


00263 

244 


00264 

244 

20S 

00266 

244 

203 

00270 

244 


00273 

244 


00274 

244 

206 


VCC <NBEST)=VEC <NBEST)-1 EXCEPT EOR LAST COMPONENT 
ALCO«nwM fOR ^eXT UMiaUE INTECER COBINATION 

DO 2D J=1,NBEST 6IVE HE NEXT UNKJUE CCWBINATICN 
42=hBEST~J41 tNCKEMENT ^BEST-TH CCNPCNENT HOST CfTEN 
JS=J24l TELLS »«: IF I INCREMENTED T«E ^©EST*TM CCHFCNENT LAST 
VECiJ2)=VECfJ2)*l INCREMEWT JZ*TH CCMFONENT 
IFCJS.6T*NBEST) CO TO 40 is THIS ^eE5T•TH CCHFONENT 
DO SO JJ=JS,N8EST NO, INITIALIZE NEXT LEVEL 

VECUJ)=VEC(JJ-1>*1 SET HIGHER ORDER TO LCLER CRCER COMP. +1 
IFIVECIhBEST) .LE.NMI CO TO 50 IS THIS LEGAL ca«0. 

CCNTINUE NO, INCREMENT THE NEXT LCtER CCMrCfJENT . 

RETURN UNIQUE CCMDINATICNS HAVE BEEN EXHAUSTED 
CONTINUE CONNECTOR, 1 HAVE A iWiaLC SLB3ET 
INCLUDE INVR.LIST MATRIX INVERSION CCCE 

OOFLfTES INVERSE FOR ‘NCL* MATRICES BV OCKDERING TECHNIQLE 

THIS ALGCRITHM PERFORMS MATRIX INVERSION CN A STPWETRIC MATRIX. 


TVE INVERSE OF "ZCOV" IS STORED IN'ZINV IN LCLCH TRIANGLE FORM. 

THIS FfiCCECURE RESTRICTS ITSELF TO T>E: ZOO SERIES 
AI^ TO INTERNAL FLAGS STARTING WITH "X* 

1 DO 221 XCL=1,NCL reduce COV. MATRIX FOR EACH CLASS 
XCT=0 COUNTS VHERE I AM IN REDUCED COV. MATRIX 
DO Z21 XC=1,NBCST REDUCED MATRIX CF RAW •^BEST‘ 

XIsVECIxC) RETRIEVE LOGICAL R^^ IM5EX 

DO 221 XD=1,XC COLUPN Ih€)EX FOR REDUCED LCVER TRIANGLE MATRIX 
XJ=VEC<XD) RETRIEVE LOGICAL CCLU44 IKCEX 

j^lly)Q^j(l^(xI~l>/24xJ CONVERT TO LO'ER TRIANGLE IKCEXING SCHE^E 
XCT-XCT^l BLM= POINTER IN REDUCED CCWARIANCE ARRAY 

ZINV<XCT,XCU=2COV<XlfC,XCU REDUCE COVARIANCE MATRIX FOR EACH CLASS 
W OOFUTE INVERSE 

DO 222 XCL=1,NCL COMMUTE INVERSE FOR EACH CLASS 
XF=1 FLAGS FIRST ROW ERROR 

IFtABStZlNVIl.xCUI.LT. .OOOOl) CO TO 230 WILL 1 BE 
VICING BY ZERO 

7INV«l,XCL)=1.0/21NVa,XCL) NO, TAKE RECIPROCAL 
DO 201 XK=2,r®EST IM3EX THRU RCV5 
DO 202 XI=1,XK INITIALIZE W>«K ARRAY TO ZERO 
zvtxi,xcu)=o 

XJ4=xk~ 1 IPPER RCW LIMIT 
OO 203 Xl=l,XM IW3EX THRU RCVS 

DO 204 XJ=1,XI IlCex THRO O0HAt« 

XIJ=(XI-l)4XI/2+XJ ODHPUTE LCU5R TRIANGLE ITCEX FOR (I.J) 

JCKJ=txK-l)4XK/2+XJ DITTO FOR (K,J) 

ZVtXl ,XCL)=2 V(XI ,XCL)*ZINV(x1J,XCL)4ZINV(XKJ,XCL) 

IF<XI-x;-:i 203,203,205 AT-J I TV.TU OCirUTINO LCRX 
XL-XI+1 »3, COMPUTE LOWER LIMIT AND CCNTINl^ 

DO 20a XJ=XL.XM MULTIPLY THE XK-TH t5CW BY THE XJ'TH COLUMN 
XKJ=(XK-1)*XK/2>XJ COMPUTE LOWER TRIANGLE INDEX FOR {K,4> 

XJl=(XJ-I>4XJ/2*Xl DITTO FOR <4, I) 

ZVIXI , XCL) =2V (XI , XCL) tZINV (XK4 ,XCL J «INV (X4 1 ,XCL) 

CONTINUE CO TO NEXT RCW 

DO 206 X4=1,XH CO THRU X4*TH COLUMN OF THE XK'TH ROW 
XK4=<XK-l)4XK/2*X4 CCWPUTE LOWER TRIANGLE INDEX (LTD FOR (K,4) 

IV (XK , XCU = ZV tXK , XCL > *2 1 NV ( XK 4 . XCL ) *Z V ( X 4 , XCL > 


COLUMN WORK VECTOR 
’.ECTCR 


COLUMN WORK VECTOR 


XK-TH ROW CF COLUMN VECTR 



Q02T4 

24* 


00277 

24* 


00300 

24* 


00300 

24* 

C BE 

00302 

244 


00303 

24* 


00306 

24* 


00311 

24* 


00312 

24* 


00313 

24* 

207 

00316 

24* 


00321 

24* 


00322 

24* 

209 

00324 

24* 

201 

00326 

24* 

222 

00330 

24* 


00330 

24* 

C 

00330 

24* 

C Ffi 

00930 

24* 

C 

00331 

24* 

230 

00333 

24* 

231 

00333 

24* 


00334 

24* 


00336 

24* 

232 

□0336 

24* 


00337 

24* 


OQ34S 

24* 

233 


CX»4C 

00347 

0035 £ 

C05»3 

00356 

oossr 

00360 

00361 
D03«7 

00362 
0O3C4 

00372 

00373 
00375 

00400 

00401 

00404 

00405 
00410 

00410 

00411 
00414 

00414 

00415 
00415 

00423 

00424 

00425 

00433 

00434 


(xKi^l)4XK/2 LTl FOR <K,M 

XF=2 Fl-Ats ERROR OCCURRED ON OTHER THAN FIRST RCVf 
lF<AB5(2INVtXRRiXCL)-ZV<XK,XCU>) *LT. .OOOOt) CO TO 230 WlU. 1 
fiC DIVIDINC 6 t zero 

ZlNV(XKR,XCU>=l.Q/tZINV(XRK,xCL)-ZV«lt,XCU) INVERSE CF EUEfCNT «,R) 

OO 207 XJ=1,XM UCex THRU R04S 
DO 207 Xl=t,XJ INDEX THRU C0LIA«S 
XJla iXJ-l)*XJ/2*Xl LTl FOR IJ|I) 

»K=tXX+l)*XK/Z LTI FCR (KiR> 

07 itNvlXJI ,XCL)=ZINV<xJl ,XCL) ♦ZINVI)0^K,XCL>*ZV(X1 ,XCL)*ZV(XJ,XCL> (J»D tNV 

DO 209 Xl=l,XM INVERSE CF THE XR-TH RCW A7C XI "TH CCLCMNS 
XKI= tXX-l)4XX/2*Xl LTI FCR (K.D 

09 21NV<XKI,XCL>=-ZlNV(XKK,XCL)#ZV{Xt,XCL) (K.I) INVERSE 

Ot CCMTIMJE CO TO ICXT RCW 
r22 OOVTINUE CO TO NEXT CLASS 
GO TO 240 SKIP ERROR PRINT 

: PRCSIXMS OH attempting to INVERT MATRIX 

rao VRtTEtPRIMT,23l) 

!31 FWMATOHl.-^W^JERSCR-W^i^ CN ATTEMPT TO COMPUTE INVERSE*, 

♦“ CF THE FOLLOWING SlB”MATRIX“ ,//) 

V«ITE (PRINT, 232) 

J32 FORMAT (• THE DIAGONAL ELEfCNTS CAm>T BE ZERO.*,/, 

TO TWO ROWS CAN BE THE SATE.*,//) 

W?17E (PRINT ,233) (VEC (XC> ,XC=1 ,^B£ST) PRINT FEATl^ES 

233 FCRMAT(1X,I6,TI15) 

Xl.=0 OOtKTS ELETFNTS IN LCv>ER TRIANGLE MATRIX 
DO 234 XC=1,M)EST IMJEX THRO NO. CHANhELS IN SLSSET 
XlaVEC(XC) GET RCW I)CEX IN BIG COV* MATRIX 
DO 255 XC=1,AC INDEX THRU CCLU+ii 
XJ=VEC<XC) CCT COL. INCEX IN BIG COV. MATRIX 

)/Z-4-XJ OCT'Fl/TE VHERE R04,CCC. ARE IN LCVCR TRIANGLE 
EUNF counter FOR TO. ELE“-ENTS IN LCVER TRIANGLE 
FI F na_l -XIU3 SAVE LdER TRIANGLE ITCEX 

235 ZPRNT(XD) =ZCOV (XINC,XCL) RETRIEVE NEXT 
C ELEMENT IN ROW FRCM BIG MATRIX 

(PRINT, 236) (ZPRNT (XM) ,XM=1 ,XC) PRINT THE RCW 

23S FCRMATUX.eElS.S) 

234 OCNTINUE CO PRINT NEXT RCW 

IF(XF.EQ.I) WRITE (PRINT, 237) TEST ERROR FLAG 

237 FORMAT (*0 ERROR OCCURRED CN FIRST RCW ,//) 

jp(xF.EQ.2) VRITE (FRINT,236) TEST ERrcCR FLAG 
23® FORMAT ("0 ERROR DID NOT OCCLR CN FIRST RCW*,//) 

VR I TE (PRINT . 239) XCL 

239 FORMAT!* THIS SUB-MATRtX IS LOCATED IN THE*,13,3H TH, 

COVARIANCE MATRIX.*) 

VftlTE (PRINT, 241 ) CHAR (XCL) 

£41 FORMAT r THIS COVARIANCE MATRIX IS FOR CLASS *,A2, 

♦- AND IS PRINTED AS IT WAS INFUT-*./) 

V«ITC<PRINT,243) (ELE(XM),XM=1,XU PRINT LCfcCR 
C TRIANGLE INDEXES 


FC*HAT(- check elements-*, 2514) 
XPTS=(NH*1)4NH/2 no. PTS IN BIG MATRIX 
write (PRINT, 244) (ZCOV{XM.XCL),XM=l,XPTS) 
FORMAT (6X , 5E1 9 . ®) 
return O FATAL ERROR. GET 0(XP. 


PRINT BIG MATRIX 
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t)04S9 

00435 

oo4se 

0D436 

00436 
00436 
00436 
00436 
00436 
00436 
00436 
00436 
00436 
00436 
00436 
OQ43T 
00440 
00443 


0O44T 

00450 

00451 
0CM92 

00459 

00460 

00461 

00462 

00463 

00464 

00469 
00467 

00470 

00471 
00471 
00471 
00474 
00479 
00477 

00901 
00501 

00902 
00910 
00511 
00511 

0091 6 
00316 

00917 
00921 
0U9E2 
00523 
00923 

00926 

00927 
00930 
00531 
00934 
00939 


24* 
246 

2S« 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

266 

RT6 

266 

206 

306 

916 

326 


240 CONTINUE CONNECTOR TO SKIP ERROR PRINT 
END 

INCLUDE DU.LIST 

C NCW SEE IF THIS 13 A COHBINATION THAT SHOIA^ BE SAVED 

C THIS PROCEDURE COMPUTES INTERCLASS OIWEKCCNCE AS FOLLOWS. ». 

C IM-ICAN VECTOR OF I’TM CLASS 

C OIJ-INTERCLASS DIVERGENCE BETVEEM CLASSSES I M3 J 
C TR-TRACE. IE. SUM OF DIAGONAL ELEMENTS 
C HC-CCVARIANCE MATRIX OF 1"TH CLASS 
C JKIN-INVERSE covariance MATRIX OF J“TM CLASS 
C N-RAMC CF MATRIX 

C OIJ=TRUK6JKlN 6 JK6IMN) -26N +TR< tIKIM*JRIN)4.<IM-JM>6<IH-JM) > 

C THIS PROCEDURE RESTRICTS ITSELF TO STATEKENTS IN THE 300 SERIES 
C AlO TO INTERNAL VARIABLES STARTING WITH W» 

WC=0 POINTS TO WHERE I AM IN INTERCLASS DIVERGENCE ARRAY 
DO 300 WA=2,NCL THE *1" CF THE OIJ CALC. IE. CLASS I. 

VAJP=uA-1 IPFER LIMIT ON LOCF 

DO 300 VB=1,WUP THE *J* CF THE OIJ CALC. IE. CLASS J. 

VC7WC61 READY FOR tCXT DU CALC. 

ZDtJ(WC>=0 INITIALIZE SCM 

yp— n VCKK CELL TO SUM IF C7E IMTERCLASS DIVERGENCE 

DO 301 WI=1,NBEST LOGICAL ROWS CF CCV. AJC INVERSE COW. MATRIX 

DO 301 WJ=t |WI LOGICAL CCLtFWS 

GET ROW I^CEX IN BIG MATRIX 

V^cVEC TWJI GET OCLIWH It^CEX 

M«.= <VK-1)6VIC/26WL ItXiEx IN COVARIANCE MATRIX CNOT REDUCED) 

WIJ=fWI-l)4WI/Z+WJ IM3EX IN INVERSE MATRIX (REDUCED) 

1 cFF DIAGONAL itKMS ARE ACC^D TWICE 
IFiWI.EO.WJ) 2=1.0 ON OIACCNAL TERMS ARE ADDED OHCE 

2j^ZD+262COV(WU-,WA)6ZINV <WI J ,WB) 7R (IK6JKIN) 

n>=ZC+26ZCOV(WU-,VvB)62INV<WI J,wA) TR<JK6lRIN) 

301 2D=ZDH-Z6(ZINV(WtJ,WA) 6ZlNVtWIJ|Vfl) ) 

16(2MEANCWt,WA)-Z>E:AN(Wt,VB)) 

26 (Z^FAN (U- ,WA) -Z7CAN (W. .WB> ) TR ( C1MN6 JKIN) 6 ( IM- JM) 6 (IM-JM) ) 

20U(WC)=ZI>-ZN«EST "2N 

IF(Z01J(WC).GT.0) CO TO 300 IS DIVERGENCE POSITIVE 
WRITE(FR1MT,310) NO, V£ HAVE ILLEGAL VALCC 

310 FORMAT (IHl , *666CRR'CR666 HAVE CCMFUTED AN ILLEGAL VALIE , 

6- FOR INTERCLASS DIVERGENCE— EITHER FCCATIVE OR ZERO.",/) 

WRITE (PRINT, 311) (VEC <WH) ,WH=1 .)eEST) PRINT CHAf-TiELS 

311 FORMAT <* USING FEATLRE StfiSET=" , 12I3,/> 

WRITE (PRINT, 512) CHAR (WA) , CHAR (VB) ,ZDU(WO PRINT CLASS 

C PAIR AM? OIJ 


312 


313 


300 

END 


.•) 


format (“ THE INTERCLASS DIVERGENCE COMPUTED fCR*, 

6” CLASS PAIR “,2A1,* IS*,E15.6,/) 
wme (PRINT, 313) ^ 

FCRMAK- ri-EASE CHECK THE STA7I37.CS . CR ,.w- -LA 

RETURN O FATAL ERROR. GET ClFP. 

CONTINIE CO COMPUTE OIJ FOR NEXT COF®INATICN CF CLASSES 

PAIRS=(NCL-l)6NCL/2 W. CF IWIQlE PAIRS 
2PCWPO INITIALIZE PACKED CHANNEL WDRO 
L=1 TtCEX THRO CHANNEL SET 

TO 62 J=1,NM INDEX THRU TOTAL »0. CF CHANNELS 
JJ=IABS(J-1) POINTS TO BIT FOSITICM IN ZPCW. INSLRE AGAINST 
tFCJ.LT.VCC(L)) CO TO 62 IS THIS CHAfWEL NOT IN THIS COMBO. 



JONES 
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oossr 

S3* 


00341 

34* 


00342 

33* 


00344 

36* 

82 

00344 

37* 

C « 

00346 

38* 

83 

00347 

39* 


00330 

40* 


00533 

41* 


00555 

42* 


00337 

43* 

no 

□0381 

44* 


00361 

45* 

C N 

00562 

46* 


00363 

47* 


00365 

48* 

C 0 

00567 

49* 


00570 

SO* 


00571 

31* 


00372 

52* 

85 

00374 

33* 


00577 

54* 

81 

00602 

S3* 


00603 

36* 


00606 

57* 

84 

00606 

38* 

C < 

00611 

59* 


00614 

60* 

ee 

00616 

61* 


0061 7 

62* 


00620 

63* 


00623 

64* 


00624 

65* 


0062T 

66* 


00630 

67* 


00631 

68* 

87 

00634 

69* 


00637 

70* 

88 

00641 

71* 


00644 

72* 


00646 

73* 


CX3647 

74* 


00650 

75* 

89 

00630 

76* 

C 

00630 

77+ 

C 

CX7632 

78* 


00654 

79* 


130657 

60* 


00661 

81* 


00662 

62* 


00663 

83* 


00666 

84* 

92 

00870 

83* 


00671 

86* 


00672 

•7* 


00673 

88* 

91 

00673 

89* 

90 


tP(j.ea.VEC(L)) FL0<iJ.t,2PCW>-t NO, URN CN BIT 
l.sL+1 LOOK AT ^CXT CHANNEL IN SET 

IP CO TO B1 HAVE I SET FUAO FOR ALL CHANNELS IN SET 

CCNTtNUE NO. CO TO NEXT CHANNEL 
3HPUTC CAPPED AVERAGE DIVERGENCE ' 

ZAVG=0 initialize AVO. DIVERGENCE 
ZM=K>IJ GET MAX. OIJ ALLOWED IN AVC. COMPUTATION 
DO 110 J=1 .PAIRS INDEX THRU CLASS PAIRS 

JPI2D1J <J) .LT. ZMI ZAVG=ZAVG+ZDIJ (J> OIJ »LT. MAX 
IFIZDIJ <J) .GE. ZM) ZAVC=ZAVG*ZM DU .GE. TO MAX DU 
ODNTINUE CO ADC NEXT CIVEF.CEMCE 
2AVC=ZAVC/PA1RS COMPUTE AVG. DU 
□W SEE IF THIS CCpeO- YIELDS THE LARGEST DU YET FOU® 

DO 65 J=l .PAIRS IMJEX THRU THE POSSIBLE DU’S 
IFIZCIJU) .LE.ZM3U<J,T>> GO TO 85 IS THIS BIGGEST 
lU I-VE Four® 

ZK3IJU.1)=ZDIJ<J> yes, save IT 

2K5IJ(J,2)=ZAVG SAVE ASSOCIATED AVERAGE DIVERGENCE 
ZK3IJ <J.S)=ZFCW SAVE CHAM^Il. CC^CO. YIELDING THIS 
CCMTINUE TRY NEXT DU 

DO 81 J=1 .PAIRS ItX^X TVKU CLASS PAIRS 

IFC2UU(JI.GE.M5U) 2DU(J)=M3U MAX- DU ALLO-EC IN AVG. CCPFUTAtICN 
2 MIN=tOOOa. INITIALIZE TO LARGE NO. IN SEARCH FOR MlN. DU 
DO 8A J=1 .PAIRS IICEX THRU CLASS FAIRS 
IF<Z0UU>.LT.ZMIN) ZMlr*:ZDUtJ) SEARQ1 FOR MINIKLM 
3.K., M3W CCPPUTE AVERAGE BY CLASS 
DO 86 J=1 ,NCU T>«U CLASSES 
ZCAVC(J)=0 initialize FOR AVG- CCH>. 

LOVPNCL-1 NO. CF “AIRS TO SLM FOR EACH CLASS 
HH=0 IM3EX THRO CU“S FOR THIS CHANvCL CCTCO. 

DO 87 1=2. NCL THE I CF THE DU 
yp =1-1 OFFER LIMIT CN J 
DO 87 J=1,UP the J of THE DU 
GO TO FC3CT OU 

2CAVC(I>=ZCAVG(I)*ZCU(m) CCFFH/TE SLM for I’TM CLASS 

2CAVC<J)=ZCAVC(JI+ZDIJ(MM) COMMUTE SW< FCR J’TH CLASS 

DO 88 J=1,NCL THRU CLASSES 

2CAVGCJ)=ZCAVG(J)/LOW CCFFOTE AVG. DU FOR CLASS 
DO 89 J=1 ,NCL I^CEX THRO CLASSES 

IFCZCAVG IJ) .LT.ZHCCI J IJ ,1 ) ) GO TO 89 

2MC0U(J|I)=2:CAVC€J) yes. SAVE IT 

2MCDU(J,E)=ZPCW ALSO SAVE COBINATlCN CF CHAW^S THAT PRCCUCED IT 
CONTINUE CO CHECK F£XT CLASS AVG. 

CHECK IF THIS COMBO. IS IN THE TCP ’M’RT 

maximum average DIVERGENCE AM5 BY KAXlMlW MINIKIM INTERCLAS^DI VERCENCE 

IFIZAVG.LE.ZMADU INPRT.l) ) CO TO 90 

CO 91 J=1,NPRT YES, RANK IT AM3 SAVE IT 
1F(2AV&.LE.2MA0IJ<J.H> C*0 to 91 
4J=J+1 YES, COPY ALL J POSITIONS TO J*1 POSITION 

DO 92 I=f-RT,JJ,'J DU‘S THAT ARE SMALLER THAN 2AVG 

2MA0IJa,l>=iMADlJ(I-l,l» lOVE SMALLER MAX. DU’S DCVN IN RANK 
ZHACIJ<1.2)=2MADUa-1.2} ALONG WITH THEIR CHANNEL VOfiCS 
2MA0U<J,*)=ZAVC INSERT NEW MAX. AVG. 01 J AT J’TH POSITION 
2MA01 J{J,Z>=2PCW SAVE CHANNELS THAT PRCCUCED IT 
eO TO 90 GO TEST FOR NEW MAX- MIN. DU 

CONTINUE CO TEST NEXT SMALLER SAVED OIJ 

IFIZMIN.LC.ZMMOIJtNPRT.D) CO TO 93 IS THIS KAX-MIN IN TOP 


IS THIS CLASS AVC. OIJ MAX. FOU® 

si CF CHAW^S THAT PRCC 

POSITIONS AS RAfKED BY 
ISS DIVE 

IS AVC. DIJ IN TOP ’NPRT’ 

: IT 

IS THIS AVG. DIJ -GT. J’fH SAVED AVG. 


-NPRT* RANK 



D06TT 

90* 


00702 

91* 


00704 

92* 


OOTOS 

93* 


00710 

94* 


00711 

93* 

99 

CX3713 

96* 


00714 

97* 


00719 

98* 


00716 

99* 

94 

00720 

too* 

93 

00722 

lOl* 


00729 

102* 


00727 

103* 


00731 

104* 


00735 

105* 

98 

00739 

106* 


OOT36 

107* 


00741 

108* 


00743 

109* 

97 

03749 

no* 


00747 

111* 


00752 

112* 


00793 

113* 


00756 

114* 

too 

□0760 

Its* 


00762 

116* 


00763 

117* 


ro7&-5 

lie* 


0076T 

119* 


00772 

120* 

lOl 

0077* 

121* 


otooo 

122# 

102 

01002 

123# 


01003 

124* 


01004 

125* 

99 

01CK36 

126* 



DO 94 J=t,NPRT ves, RANK IT ANO INSERT IT 

IP I2HXN.LC>ZKHDI J ( J I T) ) GO TO 94 tS TMIS MIN. OU .CT. J TH SAVED MIN 

JJ=J«1 YE9, COPY ALU J RANKED ELEMENTS TO J*1*ST RANK 

DO 99 1=NPRT,JJ,“1 INDEX THRU MAx-MIN DU'S THAT ARE SMALLER THAN 2MIN 

Z»«DlJ«,n=2MMDIJ(l-l,II reduce RANK CF 1*TH MIN OIJ 

2MMDlJU,e)=2MMDIja-l,ZJ *LON0 WITH ITS ASSOCIATED CHANNELS 

ZKMOl J (J,1)=ZMIN INSERT NEW MAX-MIN 

ZMMDIJfJ,2)=ZFCW ALCN& WITH ITS ASSOCIATED CHANNELS 

fiO TO 93 CO SEE IF PERCENT SEPARATION IS TOPS 

CONTINUE 60 TEST NEXT SMALLER SAVED MIN 

|F(PCTPRT.LE.O> CO TO 96 DO I SAVE MAX. PERCENTAGE DU’S 
DO 98 J=» .FAIRS YES, INDEX THRU CLASS FAIRS FOR DIJ LIMITS 
IF (ZDI J (J) «LE.O) RETURN O NEGATIVE DIJ NOT ALL^>05 

IF<20IJ(J).0T.M)IJ) 2CIJU)=f<JlJ MAX. DIJ ALLCVED IN SEPARATION CCNP. 

IF<201J(J) .LT.2.) ZCIJ<J)=Z. MIN, DIJ ALLOUEO IN SEPARATION COMP. 

CONTINUE CHECK LIMITS ON NEXT CLASS FAIR 

2PCT=1. INITIALIZE FERCENT SEPARATION RATIO 

DO 97 J=1 .PAIRS IMJEX THRU CLASS PAIRS 

IF<2PCT.GT.1.0E*35) CO TO 96 AVOID MULTIPLY OVERFLOW 

2PCT=2PCT#(2>F'CIJ<J,PCTFRU/ZDIJ(J]) CALC, RATIO OF OIJ TO SAVED POT 
IFOFCT.GE.l) GO TO 96 IS FCT SEPARATION .GT. LOVEST RAJ« SAVED 
00 99 J=1,FCTFRT YES, SEE VNERE THIS SET SHDUJ3 BE INSERTED 
3>P0T=1 , INITIALIZE PERCENTAGE RATIO 

DO JOO JJ=1 .PAIRS irCEX THRU CLASS PAIRS 

2FCT=2FCT«M^0IJ(JJ,J)/ZDI J(JJ) CALC. FERCENT SEPARATION RATIO 
IFC2PCT.CE.1.) GO TO 99 DOES THIS SET YIELD BIGGER PERCENT SEPARATION 
JJJ=J*1 YES, COPY ALL J RAIDED SETS INTO J+1 RAM< 

^-PAIRS*-! CLASS FAIRS ♦ CHAJsJnEL VERD 

DO toi I^FOTPRT, JJJi-t IM3EX TrtxU SETTS THAT HAVE SMALLER SEPARATION 
DO lot n=t.UP COPY ALL DU'S + CHAAWEL VORD 
ZFP01J(II,1)=ZFFCIJ(II,I-1) reduce RAAK CF I-TH POT SEP. 

DO 102 11=1 .FAIRS ll'Cex THRU CLASS FAIRS 
ZM=ClJ<II,J)=ZDIJ(n) INSERT rCW FERCENT SEPARATICN SET 
Z>f=OIJ<PAIRS*l,J>=ZFCVf ALONG WITH ITS ASSOCIATED CHAN^ SET 
GO TO 96 FINISHED CHECKS. CO GET ANOTHER UNIQIE: QW^EL CCMSO. 

CONTINUE CO CALC. POT SEF RATIO FOR rCxT LCl^ SET 

ETC 


OC) OF COMPILATION 


WD DIACMDSTICS. 



DISPLA 



roRtS oispla.oispl* 

rOR »*L-0«/t2'O«*4S <Z.O) 


ftlfiROUTKC OtSPi-A EJ^TSY POINT 001531 


«ORAG€ USeO- CODE (II 001345 0ATA(O> 000460 etA»« IMMMCNC2) OOOOQQ 


COMMON 

BLOCKS' 


0003 

UNITS 

000002 

0004 

MAX 

000005 

OOQS 

IW 

025416 

0006 

IMV 

005514 

0007 

YOU 

015735 


external REFEKENCES OH.OCK, NAIE) 


QOIO 

Oiv 

DDll 

ICRR6$ 

0012 

M£U$ 

0013 

NToeS 

0014 

NIOlft 

0015 


0016 

NSTOFS 

0017 

NCRR3S 


STCRACC ASSI«««NT 


(BLOCK. TYPE, RELATIVE LOCATICN, NAIC) 


OOOl 

OQOl 

coot 

ooot 

OODO 

OQOl 

OOOl 

OOOl 

OOOl 

OOOl 

ocxx> 

ooot 

OOCXl 

0000 

0001 
OOOl 

0000 

0001 
OOOl 

ooot 

OOOl 

0003 
OOOO 

0004 
0004 



000041 IL 
COOZ44 14L 
000231 16L 
000012 2L 
000107 205F 
000193 2150 
OOllSS 25L 
001203 20L 
000672 303L 
000366 3210 
000233 35F 
000313 374G 
000276 41F 
000313 44F 
000706 464C 
000734 5116 
000402 54F 
001261 6466 
001333 7056 
000337 73L 
OQOiZO 6L 
DOCOOO CARO 
I 000024 11 

OOOOOl KAXCLS 
000003 HSHO 


OOOl 

□00110 

lOL 

OOOl 

000203 

15L 

OOOl 

000261 

17L 

OOOl 

000340 

20L 

OOOO 

000115 

206F 

OOOl 

001122 

22L 

OOOl 

001163 

26L 

OOOl 

001213 

30L 

OOOl 

000731 

304L 

OOOO 

000173 

53F 

OOOl 

ce0453 

3516 

OOOl 

O014C6 

39L 

OOOl 

000560 

4116 

OOOl 

000647 

4446 

OOOl 

000715 

4716 

OOOO 

000560 

52F 

DOOl 

001051 

3506 

OOOl 

tiUl 305 

6636 

OOCl 

000415 

71L 

OOOl 

001402 

7326 

OOOl 

000521 

66L 


0005 1 024203 CHAR 

OCMW 000431 1NJF$ 

0004 000002 HEEST 

0004 OQOOOO MXHSMT 


OOOl 000167 IIL 
OOOl 000023 1526 
OOOl 000054 1706 

0000 000065 202F 

OCOO 000123 207F 
OQOl 001061 23L 

0001 000274 2676 


OOOl 

000535 

3O0L 

OOOl 

000770 

305L 

0001 

OD0431 

3366 

OOOO 

000257 

36F 

0001 

□01135 

4L 

OOOO 

000303 

42F 

OOOO 

000317 

45F 

OOOO 

000337 

48F 

OOOl 

001004 

5246 

OOOO 

000410 

56F 

OOCl 

001322 

6766 

□□01 

001344 

7116 

OOOO 

000130 

74F 

OOOl 

DOOl 62 

9U 

0005 

I 024436 

CU 

OOOO 

1 000023 

J 

0005 

024244 

PCIJ 

0003 

1 024202 

rcesT 


OOOl 

000177 

12L 

OOOO 

030057 

155F 

OOOl 

000323 

16L 

OOOO 

000073 

203F 

OOOl 

001026 

209L 

OO'.' •, 

001233 

24L 

OOOl 

001171 

27L 

OOOl 

000574 

301L 

OOOO 

OOOl 37 

31F 

OOOO 

000214 

34F 

OOOl 

0004 72 

3616 

OOOO 

000421 

40F 

□001 

000610 

4246 

OOOl 

000656 

4516 

OOOl 

000031 

5L 

OOOO 

000367 

53F 

OOOl 

001070 

5616 

OOOl 

000103 

7L 

OOOl 

001041 

72L 


OOOl 001422 7436 
DODO D00415 96F 
OOOO X D00Q31 I 
OOOO 1 000025 K 
OOOO I 000030 MM 
0003 I D242D1 NCL 


ooot D00124 13L 

0000 000033 137F 

0001 000330 19L 

0000 OOOlOl 204F 
OQOl 000265 21L 

0001 000213 2426 

OOOl 001177 20L 
OOOl 000633 302L 

0000 000154 32F 

0001 000444 3456 

OOOl 001117 STL 
OOOl 000551 4046 
OOOl 000617 4316 
OODO 000327 46F 
OOOl 000743 5046 
0001 001013 3316 

OOOl 000045 6L 
DOOl 0O1G55 70L 
OOOl 001363 7216 
OOOl 001503 7716 
OOOO 000417 99F 
OOOO 1 D00022 II 
OOOO 1 000027 LCW 
0004 000004 KPRNT 

0003 I 024200 NH 
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0005 1 024245 M’RT 
0005 I DOOOOl PRIKT 
0007 R 015672 ZCAVC 
0007 R 003700 /MCOIJ 
0007 R 004120 2MFD14 


0005 1 024245 NSHCW 
0005 I 024246 SHOW 
0005 001400 2COV 

0007 R 000760 ZHDIJ 
OOOO R OOOODO ZNCLA 


0000 I 000026 P 
OOOO 1 000052 UP 
0007 R DQOOOQ ZDIJ 
0005 OOOQQQ ZMCAN 
0007 R 015732 ZPCT 


OOOO I 000020 PAIRS 
0006 I 005500 VEC 

0006 DQOOOQ ZINV 

0007 R 013734 ZHIN 

0006 004700 ZV 


OOOO I 000021 PR IN 
0007 R 015735 ZAVC 
0007 R 004000 ZMAOIJ 
□007 R 004050 2MMDIJ 


OOlOl 

14 

ooiot 

24 

ootot 

34 

ooiot 

44 

00105 

54 

00105 

54 

D0104 

S* 

00105 

54 

00105 

54 

00106 

5# 

00107 

54 

OOltO 

34 

OOlll 

54 

00112 

54 

00115 

54 

00114 

54 

00114 

54 

00115 

54 

00116 

54 

00116 

54 

00117 

54 

00120 

54 

00121 

54 

00122 

54 

00122 

54 

00125 

54 

00123 

54 

00124 

54 

00125 

54 

00126 

54 

D0127 

54 

00127 

54 

00130 

54 

00131 

54 

00132 

54 

00153 

54 

00154 

54 

00135 

54 

00156 

54 

00157 

54 

00140 

5* 

00141 

54 

00142 

54 

00142 

54 

00143 

64 

00144 

74 

00145 

64 

00147 

•4 


subroutine OISPLA 

JGTES 

c WIS StBROirrlNE RECCPFUTES DIVER6ENCE FOR CHANNEL COBINATICNS 

C DOERHIKeD TO BE IMFCfiTAhfT AM3 OLfTFUTS THEM TO U4|T *FRINT‘ 

INCLUCE SFEC»LI5T PARAfCTERS»OITOISICNStCCMMCN 

JOCS 

ll«^LtCIT lMTEGO?(A-Y> ALL INTEGER EXCEPT Z 
PARATCTER MAXCLA=52 MAX. NO. CF CLASSES. CAN 

C INCREASE IF I^BSED. . ^ 

PARAXerTER MXMSMA=24 MAX NO. CF MEASUREMENTS <OW2CLS) . CAN INCREASE. 
PARAKETER »©ESTA=12 MAX. MO. CF PeASLREMENTS IN ‘BEST* SCBSET 
parameter RsUMA^CMEESTA+D^FSESTA/Z SLM CF 1,2,... tNO. IN "BEST S16SET 

parameter SUMA= <HXMSMA*1)4MXM1SMA/2 SIH CF 1,2,..., MAX NO. MEASLEEMENTS 
OCMMO^UMITS/CARD, PRINT 1/0 LNITS-CARD READER, PRINTER 

CCMEN /KAX/MXHSMTT,MAXCLS.K3EST,MSHO,MERMfr MAXIMA 

CCM^ON /IMF/ ZMEAN,2COW,M34,NCL, 

*^BEST,CHAR,NFRT,^CIJ,NSHCW,SHCW,DIJ IMFUT 

OIMENSICN ZMEANiMXMSMA.MAXCLA) TCLCS MEANS FOR EACH MEAStFEMCNT. CLASS 
OIMENSICN ZCCNISUMA.MAXCLA) MOLDS CCNARIANCES (LCXfR TRIANGLE) 

^ FOR EACH CLASS. 

parameter PCTFRT = 10 mo. CF COfOO-S RAMMED BY PERCENT TO BE PRINTED 
OIMEfOlCN VEC(MCESTA) HOLDS UNIQUE CCMelNATiCN CF FEATURES (CHAMNELS) 
OOMOON /IM4V/2INV,ZV,veC INVERSICW ARRAYS 

DIMENSION ZINV(RSL«A,MAxCLA> MELDS INVERSE COVARIANCE MATRICES 

FCR EACH CLASS 

^ DIMENSION ZVIMBESTA.MAXCLA) WORKING ARRAY USED WHEN COFUTING 

INVERSE BY BCRDERIM^O METMCC. 

parameter CSUMt={M<AXCLA-l)4«AXCLA/2 *MAXCLA* CC*EO*S TAKEN TVO AT A TIME 
0IMEM4SI0N ZDIJ tDSCW) VCRKIN& 6LFFER FCR INTERCLASS DIVERGENCE 
OIMENSICN CHAR(MAXCLA) HOLDS CHARACTER REFRESENTATICN CF CLASSES 
O0M04CN /YDU/ ZCl J.ZMJU rZM^CDI J ,2MACI J , 

♦ZMM3IJ,2MrCIJ,ZCAV&,ZFCT.2AVG,ZMlN FOR DU 

parameter M^SB0^^1O MWXIMLM NO. CF SHZW REOLESTS MCNORED 
OIMENSICN SHOW (MBESTA .MtSHCW) HOLDS CHANNEL CCMEO*S REQLESTEB 
parameter MFRINT=20 maximum M0. CF PRINT REQIESTS TO BE HCWORED 
DIMENSION ZMOU<CSU4,3) MCLCS MAX. CU.AVG. DIJ, PACKED CHAMWEL VCRO(PCW) 
dimension ZMCDU <MAXCLA,2> holds MWX. DU FOR CLASS )PCW 

DIMENSICX 2MADIJ <M=RINT,2) MO-DS MAX. AVG, DIJ.PCW FOR MO. PRINT REQLEST 

dimension ZnMOIJlMPRINT.2) MELDS MWX. HIN. DiJ.rCW FCH MO. PRINT REQlEST 

parameter CSUM1=CSLW+1 CLASSES TAKEN 2 AT A TIME ♦ I W>ID 
DIMENSIOJ 7MFDU<DSUM1.FCTFRTM du-s ranked by percent, pcw 
dimension ZCAVGIMAXCLA) WCKK array used to FIMC LARGEST CLASS AVG. DU 
dimension OUTDSUM) MELOS CHARACTER COMBO'S FOR CU PAIRS 


OIMENSICN ZNCLA06) HOLDS MIAXIMOM AvG. INTERCLASS DIVERGENCE 
PAIRS=TNCL“I>**^^''^2 CCMPUTE mo. OF CLASS PAIRS 
IFINSHCU.EQ.O) CO TO 1 ANY SHCW CARDS 
Pr,N= 1 YeS,FUAC telling to print 'SHOT COHBINATICNS 
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001 9Q 
00191 
W1S4 
0O15C 
ODIST 
00160 
00161 

00165 
00164 
00169 

00166 
00167 
00173 
00179 
00175 
DCIITT 
00200 
00201 

00209 

00204 

00205 

00206 

00210 
ooeli 
00212 

00219 
00214 
00217 

00220 
00222 

00224 

00225 

00226 

00290 

00291 

00292 
00299 

00295 

00296 

00297 

00240 

00241 
00244 
00249 
00247 

00251 

00252 
00259 

00259 
S02SC 
00257 

00260 
00262 
00269 
00264 
00269 
D0266 


104 

11 * 

12* 

19 * 

14* 

19 * 

16* 

IT* 

IS* 

lO* 

20* 

21* 

22* 

29* 

24 * 

25 * 

26 * 

27 * 

26* 

29* 

90 * 

91* 

52* 

93* 

94 * 

95* 

96 * 

37* 

a** 

99* 

40* 

41# 

42* 

43* 

45* 

46* 

47* 

46 * 

49* 

90* 

51* 

92 * 

99 * 

94* 

95* 

96* 

57* 

SB* 

99 * 

60* 

Cl* 

62 * 

69* 

64* 

69* 

66 * 


U=l POINTS AT -SHCW- CARO 1 

OO 9 J=1,N8EST INDEX THRU ‘SHOW" CHAI««;LS 

VCCI.»=SHOW(J.tn G€T CHANNELS FOR 01 V. COMPUTATICN 

f*i I DXV COMPUTE DlVERCENCE 

CO TO 4 PRINT OUT RESULTS 

U=tl4l BU1P SHOW CARD COUNTER 

IFai.6T.NSHCW) CO TO 1 AM I THRU WITH •SMCVT CARDS 

CO TO 2 NO. CET NEXT ONE 

11*1 COUNTS MAXIMUM INTERCLASS DIV. AVC. 

PR1N=2 FLAt TELLS TO PRINT MAX. INTERCLASS OIV. AVG. 

IJ*l COUNTS CHAr^JELS IN SUBSET FOUM3 
DO 7 J=1 ,NM II«EX THRU TOTAL MO. CHAf'ffCLS 
K*IABS<J-1) AVOID ^EGATIVE 2ERO 

IF<rLD(R,1.2MADIJCII,2)) .Ea.O) GO TO 7 I5 THIS MAX. CIV. AVG. CHA^^^EL 
IFdJ.GT.fBEST) RETURN D YES, ERROR IF MORE CHANNELS THAN IN SUBSET 
VCC(IJ)*J SAVE CHAN^ CF INTEREST FOR DIVERGENCE CALC. 

IJ=IJ41 BUM= CMANhEL COUNTER 

CONTINUE CO SEE IF hCxT CHAT2E3. IN SUBSET 

CALL DIV COMFITTE DIVERGENCE 

CO TO 4 PRINT RESULTS 

|I=ll4l BLM= MAX, DU AVG. COUNTER 

IFIII.CT.IFRT) CO TO 8 AM I THRU WITH CtJ AVG. 

CO TO 6 rO.CET NEXT ONE 

PRINTS FLAG TELLS TO PRINT MAX. MIN. DU 
11*1 COUNTS MAX. MIN. DU'S PRINTED 
U=1 COUNTS CHAWELS IN SUBSET FOLM3 
DO 9 J=1,NM INDEX TT-RU TOTAL ND. CF CHAM>CLS 
K=1ABS<J-1> AVOID -O 

IFIFLOlK.l.TWnuai -F«.0) CO TO 9 IS THIS MAX. MIN. CHATfCL 
IF<U.CT.M3EST> return a YES. ERRCR IF MORE CHATTELS THAN IN SUBSET 
VEC€U»=J SAVE CMAf2CLS CF INTEREST FCR DIV. CALC. 

1J=U41 BUMP CHAMsEL COUNTER 
CONTINUE OO SEE IF NEXT CMA^^EL IN SUBSET 
f 1 I DIV OOMPl/TE DIVERGENCE 
60 TO 4 PRINT RESULTS 

11=11+1 BUM= max. MIN. PRINT CCDUNTER 

fpilf .6T.WRT) CO TO 12 AM I THRU WITH MAX. MIN. 

CO TO 13 NO, CO GET KEXT OE 

PR1N=4 FLAG TELLS HE I AM ON MAX. PERCENT 

11=1 ooums WHICH MAX. PERCENT CN 

1J=1 COUNTS CHANTS IN SUBSET FCXPC 

DO 14 J=1,NM IM>EX THRU TOTAL NO. CF CHANNELS 

K=lAeS<J-l) AVOID -O 

lF(FLD(K,l,ZMrDU<PAlKS+l,Iin -EO-O) GO TO 14 IS THIS MAX. PERCENT CH. 
IFdJ.CT.NSEsn RETURN 0 YES. ERROR IF ICRE CHANNELS THAN IN SIBSET 
VECfIJ>=J SAVE CHANNELS CF INTEREST FOR DIV. CALC. 

IJ=U+1 BUWP CHAMVEL COUNTER 

CONTINUE CO SEE IF NEXT CHAM4EL IN SlBSET 
CALL DIV CC+PUTE DIVERGENCE 
CO TO 4 PRINT RESULTS 

I1 = XI*1 BUMP MAX. FERCE>rr PRINT CCXmER 

IFIII.GT.NFRT) CO TO 17 AM I THRU WITH MAX. PERCENT PRINT 

CO TO 19 NO, CO CET NEXT CNE 

XI=1 COUNTS WHICH CLASS AVERAGE PRINT IN CN 

TELLS HE TO PkINT CIJ AVERAGES BY CLASS 
IJsl COUNTS CHANNELS IN SUBSET POUND 
DO le J=1,NM INDEX THRU TOTAL MO. OF CHANNELS 



00?Tt 

QOtTZ 

tXi2T4 

oozre 

oozrr 

OOSOQ 

0Q902 

00909 

00904 

Q030S 

0090T 

003t0 

OD312 

0091 4 

0C915 

00316 

00917 

00920 

00929 

00924 

00924 

00926 

00926 

□0990 

00391 

00332 

00394 

00335 

00340 

00342 

00349 

00344 

00947 

ocaso 

00953 

00394 

D09S9 

00360 

00369 

00965 

00367 

D0971 

00971 

00372 

00400 

00400 

00401 

00402 

00415 

00416 

00420 

00421 

00422 
00495 
OCM96 

00440 

00441 


67 * 

66* 

69* 

70 * 

71 * 

72* 

73 * 

74 * 

75* 

76 * 

77 * 

78 * 

79 * 

60 * 

61* 

82* 

83 * 

84 * 

es* 

86 * 

87 * 

86 * 

69* 

90* 

91* 

92* 

93 * 

94* 

95* 

96* 

97* 

98* 

99* 

lOO* 

lOl* 

102 * 

103* 

104* 

105* 

106* 

107 * 

108 * 

109* 

no* 

111 * 

112 * 

119* 

114* 

115 * 

116* 

117* 

118 * 

119* 

120 # 

121 * 

122* 

129* 


R=lABStJ-l> AVOID -0 

IFlFUMR.t.ZMCDlJtllf?)) .E9.0> CO TO 18 IS THIS MAX. CLASS CHANNEL 
IFCIJ.CT.I^CST) return O yes. error if more CHAMCLS than in Sl»SET 


19 


20 


73 


VECC1J)=J 

IJ=IJ*1 

ccmtinue 

CALL DIV 
CO TO 4 
IlcII+t 


79 


SAVE CHANNEL CF INTEREST FOR OIV. CALC. 

BUMP cham«:l counter 
CO SEE IF »«:xT CHANNEL IS IN SLCSET 
cc»6*ute divergence 

PRINT RESULTS 

BUMP CLASS average PRINT COtNTER 
IFtll.CT.NCU CO TO 20 AM I THRU WITH MAX. CLASS OIJ PRINT 
CO TO 21 NO, GO GET tCXT CNE 

IF(NCL.CT.16) GO TO 70 TEST MAXlKJH NO. CF CLASSES 
IF<l«EST.&T.6> CO TO 70 TEST MAXIHLH NO. CF BEST CHAhTCLS 
XI =1 COUNTS Lire: cr print on 

FRII^S FLJkGS TELLS TO PRINT MAX. AVC. INTERCLASS DIVERGENCE 
psl INITIALIZE P 

1J=1 COUNTS CHAMPS IN SUBSET FCLN3 

DO 71 J=1,NM IhCex TVRU TOTAL NO. CF CHAh»EL5 

K=IABS(J-1) AVOID -O 

1P(FVC<K,1 .ZMCDIJ (1 1 .2) ) .Ea.Ol GO TO 71 IS THIS MAX. CLASS 
CHANhO. 

IFdJ.OT.fCEST) RETURN O YES.ERiKR IF MSRE CMAMvELS 
THAN IN SUBSET 

VECUJI=J SAVE CHAN'CLS CF INTEREST FOR DIV. CAL. 

1J=IJ*1 BUM= CHAM^EL COUNTER 
CONTINUE GO SEE IF NEXT IS IN SUBSET 
CALL DIV COfFUTTE DIVERGENCE 
DO 79 J=1,NCU IMSEX THRO CLASSES 
ZCAVCU)=0 initialize class AVC. MATRIX 
LOW=NCL-l NO. CF FAIRS TO SUM FCR EACH CLASS 
t**=0 1FC5EX THRU DU'S FCR THIS CHAT7CL C0©0 
DO 80 I=2,NCL 1W5EX THRU 1‘S CF CIJ 
UP=I-1 UPPER LIMIT CN J 
DO 80 J=1,UP THE J CF THE DIJ 
GO TO FFXT DIJ 

ZCAVGa>=ZCAVC(I)+20I J (MM) CCFFUTE SUH FOR I^TH CLASS 
ZCAVO(J)=ZCAVO(J)+ZCIJ(FW) COFFUTE SUW for J*TH CLASS 
DO 81 J=l ,NCL TWU CLASSES 

ZCAVG(J)=ZCAVG(J)/LCW COFUTE AVC. CIJ FCR CLASS 
IF(lI.NE.l) CO TO 86 HAVE I PRINTED HEADING 

VRITE (PRINT, 157) ^ 

FORMAT (IHl, -MATRIX CF DIVERGENCES FCR CHAMCL SET YIELDING . 

♦" THE MAXIMUM SEPARATION CF A CLASS FRCM ALL OTHER CLASSES*) 
VRITE(PRINT,155> (CHAR (J) , J=1 ,NCL) 

FORMAT <*0 FEATURES-, 16 (5X, AD) 

C PRINT OUT BEST CHAWELS, CLASS Af® THE MAXI MW AVERAGE DIVERGENCE 
86 OO JO <300, 301, 302, 303, 304, 305) ,NDEST 

LRITEi'FSINT.EOE) (VEC iJ) ,J-1 .NSE3T) ,CHAR(F) . (ZCAVG(J> , J^i ,NCL) 
FCRMATUSX, (Ix, I2> ,1 X, A1 , 16 (IX.F5 .0) ) 

1F<P-EG.NCL> CO TO 209 
p=p*l BUMPS P 
GO TO 209 

WRITE (PRINT ,203) (VEC (J) ,J=1 ,NBEST> ,CHAR(F) , (ZCAVG(J) ,J=1,NCL) 
FORMAT <12X,2 (IX, 12) , Ix, A1 , 16 (1X,F5,0) ) 

IF(P.Ea.NCL) CO TO 209 
p=p*l BUMPS P 
GO TO 209 


80 


81 


157 


155 


SCO 

202 


301 

203 
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00442 

1246 

302 

Q04S9 

1254 

204 

(XM36 

1264 


00460 

1274 


00461 

1264 


(XV«62 

1294 

303 

00479 

1304 

209 

(XMT6 

1314 


oosoo 

1324 


00901 

1334 


00502 

1344 

304 

00919 

1354 

206 

00916 

1364 


C0920 

1374 


00521 

1384 


00922 

1394 

309 

1X1939 

1404 

207 

00936 

1414 


00940 

1424 


00941 

1434 

209 

00942 

1444 


00943 

1454 


00949 

1464 


00946 

1474 

72 

0DS54 

1484 

74 

00559 

1494 

70 

00956 

1504 


0D557 

1914 

23 

00560 

1524 


00563 

1534 


OU564 

1944 


00966 

1554 


00570 

1964 


00571 

1574 


t»972 

1584 

37 

00974 

1394 


00579 

1604 

22 

00976 

1614 


00600 

1624 


(xieai 

1634 

4 

00603 

1644 


00604 

1654 

25 

00606 

1664 

31 

00606 

167* 


00607 

1684 


C0610 

1694 

26 

00612 

1704 

32 

00612 

171# 


00613 

1724 


00614 

1734 

27 

00616 

1744 

33 

D0616 

175# 


0061 T 

1764 


00620 

1774 

28 

00622 

1704 

34 

0O62Z 

1794 


OD623 

1804 



V«ITE(PRINT,20<) (VCC(J) , J=1,NBEST) ,CHAR<P> , CZCAVC IJ) ,J=1,NCU 
FORMAT (9X,3(tX, 12} , IX , A1 , 16 (IX ,F5,Q) ) 

IFlP.ea.NCL) 60 TO 209 
psp-t-l BIHPS P 
CO TO 209 

WRITT: (PRINT .ZOS) (VEC(J) ,J=1,NBE:ST) ,CHAR(P> , (2CAV6<JJ , J=t,NCU 
FORMAT (6X.4 (lx, 12) ,1X,AI ,16(1X,F5.0)) 

IF(P.EQ.NCL) 60 TO 209 
P=P+1 BUM^S P 
CO TO 209 

WRITE (PRINT. 206) (VEC (J) , J=1 ,f®EST) .CHAR (P) , (2CAV6 (J) , J=l ,NCU 
FORMAT (3X.5 (1X,12) ,1X,A1 , 16(1X,F5.0) ) 

IF(P^EQ.NCL) CO TO 209 
P=P+1 BUMPS P 
60 TO 209 

WRITE (PRINT ,207) (VEC ( J> , J=1 ,1«EST) ,CMAR (P) , (2CAVC ( J) , J=1 ,NCL) 
FORMAT (6 (IX, 12) ,1X,A1 , 16 (1X,F5.0) ) 

IFlP.Ea.NCL) CO TO 209 
P=P+1 BUM=S P 

ZNCLA(Il)=ZCAVG(It) SAVES MAX. AVC. DU 

11=11+1 BUMPS Lire ccwter 

IF(1I,CT.NCL> GO TO 72 AM I THRU WITH MAX. AVG. OIJ PRINT 
GO TO 73 

WRITE(FRINT,T4) tZNCLAU) , J=1 ,NCU 
FCRMATC diagonal • ,4X, 16(lX,F5-0) ) 

11=1 COUNTS WHICH MAX. DIJ PRINT CN 

PRIN=6 FLAG TELLS TO PRINT MAX. DU 

U=1 COUNTS CHANNELS IN SUBSET FOLfrC 

DO 37 J=1 ,NM IMSEX THRU TOTAL M>CGR OF CHAtl<0-S 

K=IABS(J-1) AVOID -O 

lF(FLD(A,i,2rt3UtlI,3)).tU.tJ) OO TO 37 IS THIS MAX. OU CHAWeL 
IFdJ.GT.reEST) RETURN O VES, HAVE I FOUC MORE THAN POSSIBLE 
VEC(U)=J ro, SAVE CHA^^^EL IM3ICATCR 
IJ=1J+1 BUFF COUNTER FOR CHAhteL FITC 
CONTIWE GO UTPAOC t>exT CHAANO. 

GO TO 4 PRINT RESULTS 

11=11+1 BUM= MAX DIJ PRINT COUNTER 

IF (I I. GT. PA IRS) STOP AM I FINISHED 

OO TO 23 MO, CO PRINT NEXT MAX. DU 

IFdl.NE.l) CO TO 24 HAVE 1 ALREADY PRINTED HEADING 

60 TO (25,26,27,28,29,30) ,FRIN NO, GO PRINT FRCFER HDG 

WRITE (PRINT, 31 > SHCV/ CARDS 

FORMAT (IHl,- DISPLAY OF CHANNEL CCPBINATIONS REOUESTED BY 4SHCW):- , 
CARDS. *,///> 

GO TO 24 PRINT OUTT CHANNEL DATA 
WRITE (PRINT, 32) AVC. PAIRVJISE DIV. 

FORMAT(lHl,* DISPLAY OF CHANNELS RANKED ACCORDING- TO AVERAGE", 

P* ! P-WI 5E DI VEP. (PENCE . " ; ///) 

CO TO 24 PRINT CHANNEL DATA 
WRITE (PRINT, 33) MAX. MIN DU 

FCRMATdHl," display CF CHANNEL CCHBINATIONS RANKED ACCORDING", 

*• TO MINIMUM pair-wise DIVERGENCE",///) 

CO TO 24 PRINT CHAMHEL DATA 

write (PRINT, 34) MAX. PERCENT SEPARATICW 

FORMAT OHl," DISPLAY CF CHANNEL COMBINATIONS RANKED ACCORDING", 

♦- TO MAXIMUM PERCENT SEPARATION.",///) 

60 TO 24 PRINT CHANNEL DATA 
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BOC24 

tai* 

29 

0O«2G 

162« 

39 

00625 

1636 


(X1627 

164* 


00630 

165* 

30 

00632 

166* 

36 

00632 

167* 


00633 

tea* 


00635 

189* 


00637 

190* 

41 

00640 

191* 

42 

00641 

192* 


00642 

193* 

24 

00644 

194* 

36 

D06S2 

195* 

44 

00693 

196* 


00656 

197* 

45 

00697 

198* 


00662 

199* 

46 

00663 

200* 


00664 

201* 


00667 

202* 

47 

00671 

203* 


00674 

204* 

46 

00674 

205* 


00675 

206* 


00700 

HOT* 

49 

00702 

200* 


00703 

209* 


00704 

210* 


00707 

211* 


00710 

212* 


00713 

213* 


fX}714 

214* 


00719 

215* 

50 

00720 

216* 


00723 

217* 

SI 

00729 

218* 


00727 

219* 

52 

00730 

220* 


00737 

221* 

S3 

00740 

222* 


00742 

223* 

54 

00743 

224* 


00752 

229* 

56 

00753 

226* 


00796 

227* 

98 

00737 

2£G* 


00762 

229* 

99 

00763 

230* 


00763 

231* 

C 

00764 

232* 

39 

00773 

233* 

40 

00776 

294* 


00777 

. 235* 



VR1TC»PRINT,35» MAX. AVO. DU BY Ct-ASS 

FORMAT (IHl,* DtSPCAY OF CHANNEL COMBINATICHS RANKED ACCORDING*, 

TO MAXIMUM AVERAGE DIVERGENCE BY CLASS',///! 

60 TO 24 PRINT CHANNEL DATA 
UmEtPRlNT,56) MAX. DIJ 

FORMAT (IHl,* DISPLAY CF CHAItfCL CO-*:nATTONS YIELDING MAXlMlH*. 

** PAIR-Wtse DIVERGENCE <■,//> 

VRlTEtFRINT.Al) PRINT FIRST LINE OF SlCHEADING 
miTETFRlNT,42) PRINT SECCfJD LIKE CF SLBHEACING 

FORMAT!* CLASS MAX AVG *) 

format €' PAIR 01 J DU CKAI*<CLS*,/) 

60 TO 24 PRINT CHAWCL DATA 
IFtFRIN.GT.5) CO TO 39 IS THIS MAX DIJ PRINT 
VRITE(FRINT,44) tvEC(J) ,J=^,^eEST) MO, PRINT CHAWELS 
pCfiMAT !“ CHANNELS^* ,2413 , /) 

WRITE {PRINT, 43) ZAVG PRINT AVO. DIJ 
fVQf^HATI* AVERAGE PAlR-WlSE D1 VERC£NCE=* ,FS »I ,/! 

W?ITE(PRINT,46) 2MIN PRINT MIN. DIJ 

FORMAT !' MINIMLM PAIR-WISE DIVERGENCES * ,F8.1 ,/) 

ZPCT=t initialize RATIO NJ-TIFLIER TO CNE 

DO 47 J=1 ,PA1RS IM2EX THRU CLASS PAIRS 

ZPCT=ZFCT*ZDIJ(Jl/ZH=CIJtJ.t) CALC. RATIO THIS DIJ TO MAX. SET 

JJIITE (PRINT, 48) 2FCT PRINT PERCENT OF MAX 

FORMAT (' RATIO OF THIS CHAMEL SET WITH CHAf^WEL SET*, 

♦- YIELDING HAXIMU4 PERCENT S£PARATlCN=* ,F6.2,/> 

DO 49 J=i,NCL IMJEX THRU CLASSES 
ZCAVO<J>=0 INITIALIZE CLASS AVG, MATRIX 
LCW^NCL-1 MO. CF PAIRS TO SUM FOR EACH CLASS 
M4=0 irCEX THRU CIJ*S FCR THIS CMAM«L CCPCO. 

OO 50 I=2,NCL IM3EX THRU 1*S CF DIJ 
IP =1-1 UPPER LIMIT CN J 
DO 30 J=1,UP THE J CF TPC DIJ 
M4=MH>1 CO TO PCXT 01 J 

ZCAVG(I)=ZCAV&(I)+ZCIJ(PW) C»fH/TE SIW FCR I'TH CLASS 

ZCAV&(J)=ZCAVC(J)+2DIJ(HM) compute sum for J*TH CLASS 

DO 51 J=1 ,NCU THRU CLASSES 

ZCAVG(J>=ZCAVC(J)/IjCW CCPFUTE AVG. DIJ FOR CLASS 
WRITE (PRINT, 52) SUBKJG FCR CLASS AVG. PRINT 

format C* AVERAGE DIVERGENCE BY CLASS-..*) 

V«ltE (PRINT. 53) tCHAR(J).2CAVG(J).J=l,NCL) PRINT CLASS CHARACTERS. AVG. 

format (* AVERAGE INTEKCLASS DIVERGENCE FCR CLASS *,A1," = '.Fia.S) 

WRITE (PRINT, 54) 

FCRMATIIHO,” pair-wise CI VERGEPXTE*) 

WRITE (PRINT, 56) (CIJ(J) . ZCIJ (J) , J=1 .FAIRS) 

format (lx.iodx, az,"=* ,Fe.i) ) 

CO TO NEXT FACE 


SKIP A FEW Lirr:s 


ir<NCL.CT.9) VRITE(FRINT,98) 

FORMAT (1HJ,//I 
if:s‘cl.le.o) vfitTE<rr;ir.'T,20) 

FRIN finished PRINTING THIS SET. CO GET NEXT CWE 
iINT OUT CLASS PAIR, MAX DIJ, AVG DU. AND (XANNELS 

WRITE (PRINT, 40) CIJ<II),ZKDIJ(!l.n,ZMDIJ(lI,2), <VEC(J),J=1,NBEST) 
format (3X , A2,F9 . 1 ,F6 . 1 , 24 1 3 > 

60 TO 22 GO GET NEXT MAX CIJ TO PRINT 
END 


END OF COMPILATION” 


NO DIACM3STICS 



DIV 


5-7 

/ 



roR.s oiv.oiv 

fCR **L-06/12-0®’*4 tl,0> 


subroutine OlV ENTRY POINT 001261 


STORACE USED* CCtSetll 001306 DATA(0) 000500 BCAW COHMONO 000000 
COMON BLOCKS* 


0003 

INITS 

000002 

0004 

KAX 

000005 

0003 

IM> 

025416 

0006 

INV 

005514 

OOOT 

VDIJ 

015735 


external references fBLOCK, NA^el 
OOlO NlCUS 

boil Nioes 

0012 NIOIS 
OOlS |>CRR6$ 

0014 ^ERRSS 


STORAGE ASSIGNMENT (BLOCK. TTFE, RELATIVE LOCATICN, NA^E) 


OOOl 


000016 

1476 

OOOl 


000145 

2036 

OOOl 


000415 

230L 

OOOO 


000237 

237F 

OOOO 


000273 

241F 

ODDI 


000370 

2666 

OOOO 


000352 

31 2T 

OOOl 


000572 

3666 

OOOl 


000712 

4256 

DOOl 


001155 

5316 

CX3DS 


D24436 

DIJ 

CXXK3 

I 

000172 

LOW 

0004 


000004 

►FRNT 

0005 

I 

024200 

NM 

0005 


024246 

SKW 

OOOO 


000153 

V.C 

OOOO 


000164 

LKL 

OOOO 


000132 

XCL 

OOOO 


0C0144 

XU 

oooo 


000151 

XKI 


oooo X 000152 STTS 
OOOO R DOQlSr ZD 
0007 000760 ZMDIJ 

OOOT R 013732 ZPCl 


OOOl 

000030 

1536 

OOOl 

000246 

203L 

OOOO 

000176 

231 F 

oooo 

000246 

236F 

OOOO 

000312 

243F 

OOOl 

001055 

300L 

oooo 

000370 

313F 

OOOl 

000614 

3766 

OOOl 

001032 

4536 

OOOl 

001174 

5416 

DCOO 

I 000(500 

ELE 

0004 

OOOOOl 

MAXCLS 

0004 

000003 

KSHO 

0003 

024243 

hFRT 

OOOO 

1 000175 

UP 

oooo 

I 000160 

W1 

oooo 

1 000163 

V4- 

OOOO 

I 000133 

XCT 

OOGG 

I C00140 

XIKD 

OOOO 

I 000145 

XKJ 

CCOQ 

R coDies 

2 

OOOT 

R oooaao 

ZDU 

0005 

R OOOOCX3 

ZHEAN 

OOOO 

R 000116 

ZPRNT 


OOOl 

000044 

1576 

0001 

000155 

2116 

CX3O0 

000215 

232F 

0000 

000256 

239F 

oooo 

000317 

Z44F 

oooo 

000321 

31DF 

OOOl 

000452 

3176 

DOOl 

000635 

4106 

OOOl 

001075 

5026 

OOOl 

001204 

5476 

oooo 

1 000174 

I 

0004 

000002 

l«EST 

0004 

oooooo 

MXMSMT 

0005 

024245 

NSHOW 

0006 

I 005500 

VEC 

00(30 

I 000165 WU 

OOOO 

I 0C3O167 

WM 

OOOO 

1 000136 

XD 

DOOO 

t 0(30137 

XJ 

OOOO 

I 000150 

XKK 

0007 

r» OJ5733 

ZAVG 

0006 

R OOOOOO 

ZINV 


OOOT R 01ST34 2MIN 
0006 R 004700 ZV 


0001 000072 171C 

DQOl 000170 2140 

0000 000233 233F 

0001 ooozeo 240G 

DOOl 000321 253C 
0001 000443 310G 

0001 000466 323G 

OOOl 000655 414G 
OOOl 001133 5160 
0003 000000 CARD 

OOOO 000410 INJF$ 
0005 I 024244 M3IJ 
0005 I 024202 NBEST 
OOOO 1 O'lOlTQ PAIRS 
OOOO I 000154 WA 
OOOO I DOQ161 WJ 
OOOO I 000155 WJP 
OOOO 1 000141 XF 
OOOO I 000147 XJI 
OOOO X 000146 XL 
0007 R 015C72 ZCAVG 
0007 004000 2MACIJ 

0007 004050 Z»CIJ 


OOOl 000123 2000 

OOOl 000222 227G 

0000 000235 236F 

OODl 000622 240L 

0001 000336 2560 

0000 000344 31 IF 

0001 000512 335C 

OOOl 000677 422C 

OOOl 001146 5256 

0003 1 024203 CHAR 

OOOO I 000171 J 
OOOO 1 000173 MM 
0005 I 024201 NCL 
0003 I OOOOOl PRINT 
OOOO I 000156 VB 
OOOO 1 000162 WK 
OOOO 1 000134 XC 
OOOO I 000135 Xl 
OOOO 1 OOOl 42 XK 
OOOO 1 000143 XM 

p 00140^ ^COV 
0007 DQ3709 2MCDI J 

0007 R 004120 2MFD1J 



SU&ROUTt»C DIV 


JONES 


ootot 

1« 

OOlOt 

2* 

ootot 

3* 

ootot 

4« 

ootos 

S<«> 

00103 

9« 

D0104 

S* 

DD109 

9* 

ootos 

S* 

ootoe 

s* 

ootor 

9« 

ootio 

94 

CXllll 

94 

OOtl2 

94 

D01t3 

94 

00114 

94 

00114 

94 

OOltS 

94 

00116 

94 

00116 

94 

03117 

94 

OOlZO 

94 

00121 

94 

00122 

94 

00122 

94 

00123 

54 

00123 

54 

00124 

54 

00125 

94 

03126 

54 

0U127 

94 

001 27 

94 

00130 

94 

00131 

54 

00132 

54 

00133 

94 

00134 

54 

00139 

34 

00136 

54 

00137 

94 

00140 

54 

00141 

94 

00142 

94 

00142 

54 

00143 

64 

00143 

74 

00144 

64 

00144 

94 

00149 

104 

00149 

104 

00149 

104 

00149 

104 

DD14S 

104 

00149 

104 

00149 

104 

00149 

104 

D014S 

104 


C THIS Sl®KOOTlNE COMPUTES DIVERGENCE ATD PARAMETERS NEEDED ©Y DISPLAY 
C FCR A GIVEN SET CF CHANNELS 

INCLUDE SPEC .LIST PARAMETERS, DIMENSIONS. COMMON 

^ JCNES 

IhCLtClT INTEGER <A-Y) ALL INTEGER EXCEPT 2 
parameter MAXCLA=32 max. no. CP CLASSES. CAM 
C INCREASE IF NEEDED. 

PARAPETER MXMSMA=24 MAX MO. CF tCASUREMENTS tCHANNELS) . CAN INCREASE. 
PARA^ETER ►©E3TA=12 MAX. NO. CF MEASUREMENTS IN -BEST* StBSET 

PARAfCTER RSUMA= (f©ESTA*l)>W-SESTA/2 SW CF 1,2 NO. IN "BEST* SLBSET 

PARAfETER SUW= (MXMSMA+1) +MXMSMA/2 SLW CF 1,2,..., MAX NO. MEASUREMENTS 
OU<MCN/UJITS/C.*.~D .PRINT I/O LNITS'“CARC READER, PRINTER 
OCMMON /MAX/MXMSMT,MAXCLS,W3EST,MSH0,M=RNT MAXIMLWS 
COW40N /INF/ 2tEAN,2COV,NM,NCL, 

««EST,CHAR.NFRT,MDIJ.NSHCW,SHCW.DIJ IW=UT 

OITENSICN 2^EAN(MXMSMA,MAXCLA) HOLDS MEANS FOR EACH >EASU^EMENT , CLASS 
DIMENSION 2COV(SUKA,MAXCLA> HOLDS COVARIANCES (LCVER TRIANGLE) 

^ BDR EACH CLASS. 

PARAfETTER PCTFRT=1D NO. OF CCteo‘5 RAWED BY PERCENT TO BE PRINTED 
DIMENSION VECCKBESTA) HOLDS UNIOUE CCMBINATICN CF FEATURES (CHATFCLS) 
COMMON /INV/2INV,ZV,VEC INVERSION ARRAYS 

OIICNSICN ZINVtRSLMA.MAxCLA) MOLDS INVERSE COVARIANCE MATRICES 
^ FCR EACH CLASS 

DIPEHSICN ZV<^eESTA,MAXCLA) VCRKING ARRAY USED VHEN COPFUTING 

INVERSE BY BCRDERING METHCO. 

PARAFerO? DSU4=(MAxCLA-l)*«AXCLA/2 -MAxCLA* CO^S. TAKEN TVO AT A TIME 
DIKENSICN ZOIJC3SUM) WCKKING BUTER FOR INTERCLASS DIVERGENCE 
DIMENSION CHAR(HAXCUA) H3J5S CHARACTER REPRESENT AT ICN CF CLASSES 
COtlON /YDiJ/ ZDU.ZMJlJ.ZMCUU.ZMADIJ, 

#2MvCIJ,ZM=01J,ZCAV&,ZPCT,2AV&,ZMIN FCR CIJ 

PARAhETER MSHCLFIQ MAXIMUM NO. CF OXVf REOlFSTS HCNORED 
DIKENSICN SHOW t^©ESTA,MSHOW) HtXDS CHAfJCL CCPBO'S REQUESTED 
PARAJ-CTER WRINT=2a MAXIMUM NO. OF PRINT REQUESTS TO BE VCNORED 
OIIENSICN ZW?IJ<DSUM,3) HOLDS MAX. DIJ.AVG. CIJ, PACKED CHAWCL VCRDIFCW) 
ClfENSlCN ZMCDIJ<MAXCLA,2) HOLDS MAX. OIJ FOR CLASS, PCW 

OII-CNSICN 2MA0U IMPRINT, P) HOLDS MAX. AVG. CIJ, PCW FOR NO. PRINT REQUEST 

OIPENSICN ZVP0IJIM=RINT,2> HOLDS MAX. MIN. DIJ.FCW FOR MO. PRINT REQUEST 

PARAMOTER CSUM1=DSUIM+1 CLASSES TAKEN 2 AT A TIKE + 1 WORD 
DIMENSION ZKFDI J (DSUIMl .FCTPRT) DIJ*S RANKED BY PERCENT, PCW 
DIKCNSICN ZCAVG (MAXCLA) WORK ARRAY USED TO FIND LARGEST CLASS AVG. DU 
DIMENSION DIJ(CSUIM) HOLDS CH^ACTER CCPBO*S FCR DU FAIRS 

EKO 

DIMENSION ELE(RSUMA) PRINT BUFFER. HOLDS ELEMENTS 

C CF COVARIANCE MATRIX ON ERROR. 

DIMENSION ^RNTCMBESTA) PRINT BUFFER. HOLDS RCW 
f CF COV*Fl*NCE MATRIX CN FPRC*’ . 

INCLUDE INVR.LIST PROCEDURE INVERTS COVARIANCE MATRIX FOR ALL CLASSES 

JONES 

C COFUTES INVERSE FOR "NCL* MATRICES BY BORDERING TECHNIQUE 

C THIS ALGORITHM PERFORMS MATRIX INVERSION CN A SYMMETRIC MATRIX. 

C SEE 'STATISTICAL COMPUTATIONS CN A DIGITAL COMPUTER* ©Y HEMMERLE- PG 73 
C THE INPUT MATRIX IS *ZCOV* IN LCVER TRIANGLE FORM 
C THE INVERSE CF ’ZCOV* IS STORED IN-ZIHV* IN LO.ER TRIANGLE FORM. 

C THIS PROCEDURE RESTRICTS ITSELF TO THE 200 SERIES 



00145 

109 

00149 

lO* 

00151 

10* 

W1S2 

lO* 

00155 

10* 

00156 

lO* 

mi«t 

to* 

00162 

io« 

00163 

lO* 

0O194 

to* 

00194 

lO* 

1X3170 

ID* 

00173 

lO* 

00174 

io« 

00174 

io« 

00176 

io« 

00177 

lO* 

00202 

lO* 

1X3205 

10* 

DQ207 

lO* 

ooeiD 

lO* 

00213 

io« 

00219 

ID* 

00217 

ID* 

00220 

lO* 

00222 

lO* 

00229 

lO* 

00226 

ID* 

□0231 

lO* 

OOCT2 

ID* 

00233 

ID* 

00235 

ID* 

00237 

lO* 

00242 

ID* 

00243 

ID* 

00245 

ID* 

00249 

ID* 

00247 

to* 

00247 

ID* 

00251 

ID* 

00252 

ID* 

00255 

ID* 

00260 

ID* 

00261 

to* 

00262 

to* 

00265 

ID* 

00270 

ID* 

00271 

lO* 

00273 

lO* 

00275 

10* 

00277 

to* 

00277 

ID* 

00277 

lO* 

00277 

10* 

00500 

to* 

00302 

to* 

00502 

to* 


ZZl 


C Al® TO INTERNAL FLACS STARTING WITH *X* 

111 DO 221 XCL=1.NCL REDUCE COV. MATRIX fOR EACH CLASS 
XCT=Q COUNTS VAIERE 1 AM IN REDUCED COV. MATRIX 
t» 221 XC=1,N8EST REDUCED MATRIX CF RANK *r©EST* 

K1=VEC(XC) RETRIEVE LOGICAL ROW INDEX 

DO 221 XD=1,XC CCLm>« UCEx FOR REDUCED LCLER TRIANGLE MATRIX 
V j^yEC {XDI RETRIEVE LOGICAL CCLIAW IFCEX 

XlM3=Xl*tXl-l)/2+XJ CONVERT TO LCVER TRIANGLE INDEXING SCHEME 
XCT=XCT+1 eUMP POINTER IN REDUCED COVARIANCE ARRAY 

ZINV<XCT,xCL)=ZCOVtXlhO,XCL> REDUCE COVARIANCE MATRIX FOR EACH CLASS 

M3W COFUTE INVERSE „ 

DO 222 XCL=1 ,NCL COMPUTE INVERSE FCR EACH CLASS 

XFal FLAGS FIRST ROJ ERROR 

|F(ABStZINV(l.xCLl>.LT. .OOOOl) GO TO 230 WILL 1 BE 
DIVIDING OY ZERO 

ZINV<l,XCL>=1.0/2INVa,XCL) NO, TAKE RECIPROCAL 
DO 201 XK=2,ICEST IM3EX THRU RCVS 
DO 202 XI^I.XK initialize VCRK ARRAY TO ZERO 
2V€XI.XCL>=0 

XM=XK-1 OFFER RCW LIMIT 

DO 203 Xl=l ,XM IM3EX THRU ROWS 

DO 204 XJ=l,Xl IM?eX THRU CCLLFWS 

XIJ=<XI-l)#Xl/2+XJ CCMFUTE LCVfR TRIANGLE ItCEX FOR (1,J> 

*^i«UXCLrSv«I^xCL>?zlL(XlJ.XCL)«INV(xKJ,^^ CCLLFM VORK VECTOR 
205,203,205 AH I THRU COMPUTING VCRK VECTOR 
vt=xl+l ND, OCFFUTE LOWER LIMIT AJ^ CCNTlNLE 

DO 208 XJ=XL,XH MULTIPLY THE XR-TM RCW BY THE XJ*TH CCLIW^ 
j<KJ=tXK-l>*XK/2+XJ COMPUTE LCWER TRIANGLE IhCEX FOR (K.J) 
XJT=fXJ-l)#XJ/2+Xl DITTO FOR (J,l) 

ZV«XI.XCU=ZV(XI.XCL)+ZINV(XKJ,XCL>«INV<XJI,XCL> CCLW« VCRK VECTOR 

CONTINUE CO TO FEXT RCW 

DO 206 XJ=1,XM GO THRU XJ-TH CCLUT^ CF ^HE W‘TH 
«cJ-(WC-l)=frXK/2»XJ CCMFUTE LCVER TRIANGLE IICEX fl-TI) FOs (K,J> 

ZVI»,XCU=ZV<xK.XCL)*2INV(XKJ.XCL)«V<XJ,xa.) XK*TH KM CF CCLUMM VECTR 
XKK=(XK+1)*XK/2 LTI FOR <K,K) , 

jtr-2 FLAGS ERROR OCCURRED ON OTTER THAN FIRST RCW 

IFIAeStZINV(XKK.xCL)-ZV(XK.XCL».LT. .00001) GOTO 230 WILLI 

inverse ee eeeheh, <k,« 

DO 207 XJ=1,XM ITCEX THRU ROWS 

DO 207 Xl=l.XJ INDEX THRU CCLW^ 

XJI=CXJ-l)«J/2+XI LTI FOR <J.I> 

.NVERSE CF THE XX-m KW Xl'TH «XT«e 
XK1=(XK-1)*XK/2+XI LTI FOR (K,I> 

ziNVix:;; ,xcL)-r”--K:t,xcL;^2V<xl ,xcl> 

CONTINUE GO TO NEXT RCW 

continue go to next class 

CO TO 240 SKIP ERRCR PRINT 


202 


204 


203 


208 

203 


206 


207 


ED9 

201 

222 


;,1) INVERSE 


PRCBLEMS on attempting TO INVERT MATRIX 


230 

231 


14? I TE (PR INT , 231 ) 

FORMAT ATTEMPT TO CCMFUTE INVERSE*, 

CF THE FOLLOWING SL©-MATKIX* ,//> 
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00303 

ID* 


00303 

10* 

232 

00303 

IQ* 


(»30« 

lO* 


00314 

10* 

233 

00313 

to* 


00316 

ID* 


00321 

lO* 


00322 

lO* 


00323 

lO* 


00326 

lO* 


00327 

ID# 


(X3330 

to* 


00331 

io« 

235 

(X)331 

ID* 

C EL 

00333 

lO* 


00341 

lO* 

236 

00342 

lO* 

234 

00344 

lO* 


00347 

104 

237 

00350 

ID* 


00353 

ID* 

238 

00354 

to* 


00357 

lO* 

239 

00357 

to* 


DD360 

lO# 


(X»63 

104 

2.. ' 

00363 

lO* 


(XX364 

104 


CXJ364 

104 

C IF 

(X1372 

to* 

243 

00373 

104 


00374 

104 


00402 

ID* 

244 

00403 

to* 


00404 

lO* 

240 

00404 

lO* 

END 

00403 

11* 


0D40S 

11* 

C 


WRITE tFRINT,Z38) 

FORMAT f- THE DIAGONAL ELEMENTS CAWOT BE ZERO.*,/, 

*• NO TWO ROWS CAN BE THE 5Af€,“,//) 

VJiXTC IFRINT «23S) IVEC (XC> »XC=1 jNBESTl PRINT FCATW'ES 

FORMAT ax. 18, Til 5) 

)(I.=0 COUNTS ELEMENTS IN LOk£R TRIANGLE MATRIX 
DO eS4 XCil.NBEST INDEX THRU NO. CHANNELS IN SIBSET 
Xl=VEC(XC) GET ROW INDEX IN BIG COV, MATRIX 
DO 235 XD=1,XC IM3EX THRU CCUA^ 

XJsVEC(XC) GET COL. INDEX IN BIG COV. 
XlMi=XI*(XI-l)/2*XJ COMPUTE WHERE l=^'- 
XL=XL*1 COUNTER FCR NO. ELEMENT:- 

ELE<XL)sXIND SAVE LOWER TRIANGLE INDEX 
ZFRNT(XD)=2COV(Xl«®,XCL) RETRIEVE »>EXT 


rstx 

L. ARE IN LOWER TRIANGLE 
LCWER triangle 


VJmEtFRINT,236) (2FRMT <XM) ,X«=1 ,XC) PRINT THE RCW 

FORMAT ax,8E13.8) 

CONTINUE CO PRINT >CXT RCW 

IFIXF.EO.D WRITE IPRINT ,237) TEST ERROR FLAG 
FORMAT ("O CRRCR OCCURRED ON FIRST RCW” ,//) 

IF<XF.EQ.Z) Wi<ITE<FfiINT,238) TEST ERROR FLAG 
FORMAT C”0 ERRCR DID NOT OCCIR ON FIRST RCW",//) 

VR ITE (PR I NT , 239 ) XCL 

FORMAT!” THIS SUO-MATRIX IS LOCATED IN IPE* , 13,3H”TM, 
COVARIANCE MATRIX.”) 

WRITE (PRINT, HAt) CHAR(XCL) 

F(3?MAT(” THIS COVARIANCE MATRIX IS FOR CLASS ”,A2, 
AhO IS PRINTED AS IT VTAS IM^UT.”,/) 

VRITE (PRINT, 243) (ELE(XM) ,XM=1 ,XL) PRINT LCVER 


FORMAT!” CHECK ELEMENTS-” ,2514) 
jiPTS= (NM+1) >NTM/2 NO. .FTS IN BIG MATRIX 
WRITE (PRINT, 244) (2COV(XH,XCU ,XK=1,XFTS) 
FORMAT <6X,3E15.8> 

RETURN O FATAL ERRCR. GET CtPP. 

GCNTINUE CONFECTCR TO SKIP ERROR PRINT 


PRINT BIG matrix 


CX3405 

00405 

00405 

00405 

00405 

00405 

0CM05 

00405 

00403 

00405 

0O40e 

00407 

00412 

00413 
0O4J« 
00417 

00420 

00421 


11 # 

11 * 

11 * 

tl* 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 

11 * 


include OIJ.LIST FRCXrECURE TO CCFfUTE DIJ FOR ALL CLASSES 


IM-FEAN VECTCR OF I”TH CLASS 

OIJ-INTERCLASS divergence between CLAS5SE5 I AhC J 
TR-TRACE. IE. SUM CF DIAGONAL ELEMENTS 
IK-COVARIANCE MATRIX OF I”TH CLASS 

jmu — inverse covariance matrix of J”TH class 

N-RANK OF MATRIX 

DIJ=TR(1K*JKIN ♦ JK*IKIN) -2*N -►TR ( (1KIN+ JKIN) * (IH-JM) *(IM- JM) ) 
PROCEDUCE RESTRICTS ITSELF TO STATEfSINTS IN THE ”32 SERIES 
AM) TO INTERNAL VARIABLES STARTING V/ITH W. 

WC=0 POINTS TO WHERE I AM IN INTERCLASS DIVERGENCE ARRAY 
DO 300 WA=2,NCL THE ”1” CF THE DIJ CALC. IE. CLASS 1. 

VIUP=WA-1 UPPER LIMIT ON LOOP 

DO 300 VjB= 1,VA)P THE “J* OF THE DIJ CALC. IE. CLASS J. 
hC=V)C*l READY FOR NEXT CIJ CALC. 

ZDIJ(WC)sO INITIALIZE SIH 

ZD=(J VCFK CELL TO SUM UP CNE INTERCLASS DIVERGENCE 

DO sot WI=l,f«EST L(301CAL ROWS CF COV. AND INVERSE COV. MATRIX 



JONES 
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00424 

11* 

OD42T 

11* 

00430 

11* 

00451 

11* 

00432 

11* 

00433 

It* 

00454 

11* 

00456 

11* 

00437 

11* 

00440 

11* 

00440 

11* 

00440 

11* 

00443 

11* 

00444 

It* 

(XM46 

11* 

00430 

11* 

0O450 

11* 

00451 

. 11* 

00457 

11* 

004 GO 

11* 

00460 

11* 

00465 

11* 

00465 

11* 

00466 

11* 

00470 

11* 

00471 

11* 

00472 

11* 

00472 

11* 

00475 

12* 

D0475 

12* 

C0473 

12* 

00475 

12* 

00475 

12* 

00476 

12* 

00477 

12* 

00500 

12* 

00501 

12* 

00504 

12* 

00S06 

12* 

ODSIO 

12* 

00512 

12* 

00514 

12* 

00515 

12* 

00520 

12* 

00522 

12* 

00525 

12* 

00524 

12* 

00527 

12* 

00550 

12* 

00533 

12* 

00534 

12* 

00535 

12* 

00540 

12* 

00543 

12* 

D054S 

12* 

00545 

13* 

00546 

14* 



oo 301 WJ31.WI LOGICAL COLICS 
vr=VECiWl) CCT ROW INDEX IN BIG MATRIX 
LL^veCtWJ) GET COLUMN ITCEX 

MlL=fWK-lJ#VK/2^WL IMJEX IN COVARIANCE MATRIX (NOT REDUCED) 
WIJ=(W1-1>W1/2*WJ IW5EX IN INVERSE MATRIX (REDUCED) 
j_g ALL OFF OIAGCKAL TERMS ARE ADDED TWICE 
IFtWl.EQ.WJ) Z=1.Q ON DIAGONAL TERMS ARE ADDED CNCE 
2D=ZC*2*‘ZC0V CWRL,WA)*ZINV (W1 J,ve> TR(IK*JMN) 

ZD=ZD+24=ZC0V CWKL ,vB>*ZIHV{WIJ,WA) TR(JR*IKIN) 

301 2D=ZD+Z* tZlNVCWIJ,WA)+2INV(WlJ,ve)) 

!♦ (ZKEAN [W^ , WA ) -Zt^AN (LK , V« ) ) 

Z*tZMEAN(VJL»WJA>-2MEAN(WL,vB)) TR { (IMN* JMN) ♦(!«- JM)*(IM-JM) ) 

ZDl J (WC) =ZD-Z^EST -2W 

IF (ZDIJ (WC> .GT.O) CO TO 300 IS DIVERGENCE F06IT1VE 
V#TITE{FRII'fTr310) HO, WE HAVE ILLEGAL VALLE 

310 FORMAT (lHl,“’(=t=!ERRCR4=S=4‘ HAVE CCMTUTED AN ILEEGAL VALLE', 

♦■ FOR INTERCLASS DIVERGENCE — EIIFER ^EGATIvE CR ZERO.*,/) 

VRITE(FRINT,311) (VEC (UM) ,VJ<=1 ,fEEST) PRINT CHAWELS 

311 FORMATl* USING FEATURE 5LBSET=" ,1213,/) 

WRlTElFRINTiSlZ) CHAR (ViA) , CHAR (WB) ,2DIJ CWC) PRINT CLASS 

C PAIR A^D DU 

312 FCRKATC THE INTERCLASS BIVEKOENCE CCM^UTEB FOR", 

♦" CLASS FAIR *,2A1,' IS'.EIS.S,/) 

WRITE (PRINT, 513) _ 

313 FORMAT (“ PLEASE CHECK THE IM’UT STATISTICS FOR THESE CLASSES. ) 

return o fatal errcr. get CLM*. 

TOO COVTINUE CO COEUTE DIJ FCR FEfT COCINATICN CF CLASSES 


end 


include stat,list 


400 


402 


FRCCEEURE to CCMP. AVG. DIJ, MIN. BIJ, AVG. CLASS DIV 

JOES 

WIG rROCCruRE COrUTES AVERAGE FAIR-WISE DIVERGENCE, MINIMUM PAIR-WISE 
DIV., Ate AVERAGE DIV. BY CLASS. 

THIS PROCEDURE RESTRICTS ITSELF TO STATEFENTS IN WE AOO SERIES 
PAIRS=CNa.-l)<NCL/2 CCM=UTE NO. CF CLASS PAIRS 
ZAVC=0 INITIALIZE AVERAGE FAIR-WISE DIV. 

2MIN=10000Q. INITIALIZE TO LARGE MO. IN SEARCH FCK MIN. PAIR-WISE DIV. 
to 400 J=l, PAIRS irCEX THRU PAIRS 

IF(ZDU(J) .LE.D) RETURN O ERRCR IF NCNFOSITIVE CIJ 

|F(ZD1J,(J) -GT.MJU) ZCIJ(J)=rt3U MAX. DIJ ALLCLED IN AVG. CCMPUTATICN 
IF(Z0U(J>.LT.ZMINI 2MIN=ZDIJ(J) save MlNlMlM FAIR-WISE DU 
ZAVG=ZAVG+ZDIJ (J) SUM FAIR-WISE CIV. 

ZAVC=ZAVG/PA[P.S AVG. PAIR-WISE DIV, 

00401 J=1,NCL THRU CLASSES 

ZCAVG(JJ=0 initialize FCR AVG. DIV. BY CLASS CCMP. 

LOW=NCL-I WO. CF PAIRS TO SUM FCR EACH CLASS 
W4=0 IfCEX THRU FAIR-WISE CIV. 

DO 402 1=2. NCL INDEX THRU I'S CF DU 

UP=I-1 UPPER LIMIT CM J 

DO 402 J=1,UP THE J CF THE DU 
CO TO NEXT DI J 

ZCAVG(r)=2CAVCa)*2DU(MM) COMPUTE PAIR-WISE SUM FCR I"TH CLASS 

ZCAVG(J>=ZCAVG(J)*ZDU(MH) COMPUTE PAIR-WISE SlM FCR J'TH CLASS 

DO 405 J=1,NCL thru CLASSES 

ZCAVG(J)=2CAVC(J]/LCW CCHFUTE AVG. PAIR-WlSE DIV. FOR EACH CLASS 


405 

END 


ZPCT=l. INITIALIZE PERCENT SEPARATION 
OO 1 J=1 .PAIRS IlVDEX THRU CLASS PAIRS 
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DOSSt 

19« 

OOSS^ 

16* 

ODSSS 

IT* 

oossr 

16* 

00960 

19* 


irtZDlJ(J) .CT.WJtJ) ZOU<J)=MDIJ 
lFtZ01JtJ)-LT.2.) ZDIJ(J>=Z. 
IPCT=2FCT*2HP01 J { J ,1) /2DI J IJ) 


K*,:. otJ allofED in separation CCNP. 
hin cij allowed in fct separation comp. 
CALC. RATIO OF NEW CIJ TO MAX. PCT 


RETURN 

END 


CNO CF CCWlLATtCW’ 


no DiACNoeTJcs. 



Input Data Deck 

{from flight line C-1) 

his deck was input to produce the Example of Output. 
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Oata.l Oat a, 

OATAooe-Ruie«r-io os/i2-oa'4i*t2 


000001 

CHANNCLS=12 CLASSes=D9 

BEST=Q4 




000002 

ICAN 

.16997769*03 

.17495210+03 

.19327100+03 

.19252034+03 

.16906102+03 

OOQOOO 

MEAN 

.16709183403 

.19051837+03 

.17088713+03 

.18449409*03 

.17266023+03 

000004 

(CAN 

.16242322403 

.18200590403 




OOOOOS 

>CAN 

.17085855403 

.17711908+03 

.19527697+03 

.19497829+03 

.17264605+03 

oaooa« 

>EAN 

.16823223403 

.19251315+03 

.17829999+03 

,19304539+03 

.18202960+03 

000007 

»CAN 

.14818815403 

.17376513+03 




oooooe 

MEAN 

.17872218403 

.16174953403 

.19723617+03 

.19657800*03 

,17516097+03 

000009 

MEAN 

.17049099403 

.19132753+03 

.169 50404+03 

.18331013+03 

.16658421+03 

OOOOIO 

mean 

.14759664403 

.17001491+03 




OQOOtl 

»CAN 

.18260169403 

.18367144+03 

.19701955*03 

.19534094+03 

.17435737+03 

000012 

►CAN 

.17377575403 

.19161604+03 

.16234224+03 

.17112711+03 

,15899739+03 

000013 

MEAN 

.17283181403 

.18596349+03 




000014 

KriAN 

.18234438403 

.18560946+03 

.20049707+93 

.20072515+03 

.18267511+03 

000013 

►CAN 

.17498570403 

.19682326+03 

,16611500+03 

.19962033+03 

.17815334+03 

000016 

»CAN 

.11384275403 

.14711826+03 




000017 

^EAN 

.17783389403 

.18299663+03 

.19925926+03 

.19918069+03 

.17933109+03 

000018 

MEAN 

,16979460403 

.19422109+03 

.18477328+03 

.19993602+03 

.17645118+03 

000019 

ICAN 

,10342536403 

.14320426+03 




000020 

MEAN 

,17489414403 

.17520609+03 

.19134081+03 

.18964935+03 

.16111868+03 

ODD021 

MEAN 

.19829190403 

. 18242902+03 

.15211467*03 

.16892622+03 

.15519086+03 

000022 

»CAN 

.15898716403 

.17907938+03 




000025 

KEAN 

.16494736403 

.16969132+03 

.18957894*03 

.18815505+03 

.16247937+03 

000024 

lEAN 

.16552632403 

.10764068+03 

.15987767+03 

.17153627+03 

.16476387+03 

00002S 

►EAN 

.18318918403 

.19461735+03 




OOCX326 

♦CAN 

,17433800403 

.17444000+03 

,19021 799+03 

.18774299+03 

.15613600*03 

000027 

>CAN 

.15355899403 

.17782400+03 

.13898700+03 

.15473799+03 

.13909599+03 

000028 

MEAN 

.15609799403 

.17002300+03 




000029 

ODVAR 

.65274073401 

.53574336+01 

.77711774+01 

.39174090+01 

.39079335+01 

000030 

OOVAfi 

.39718175+01 

.37877160+01 

,42324261+01 

.28162173+01 

.38085137+01 

000031 

COVAR 

.74806794401 

.82846057+01 

,34437276+01 

.6tD450481+Dl 

.13914199+02 

000032 

COVAR 

.64108749+01 

.73504820+01 

.48119443+01 

.54267190+01 

,11094011+1^ 

000033 

OOVAR 

.12005847+02 

.44160372+01 

.48613533+01 

.35376967+01 

.37425675+01 

000034 

OOVAR 

.72447743+01 

.70952387+01 

.56575742+01 

.71964298+01 

.86472496+01 

000035 

COVAR 

.56510909+01 

,63311547+01 

.12734403+02 

.11478605+02 

.77670126+01 

000036 

COVAR 

.16037153+02 

.55520678+01 

.70039848+01 

.44864996+01 

.52108884+01 

00003T 

OOVAR 

.iniDCM92+02 

.87373384+01 

.62906535+01 

♦12599464+02 

.12289523+02 

000038 

OOVAR 

.57565647+01 

.70405760+01 

.48323850+01 

.53846570+01 

.1093X712+02 

000039 

OOVAR 

.10475405402 

.72931145+01 

.13130637+02 

.11170218+02 

.14280519+02 

000040 

COVAR 

.32038013401 

.23827266+01 

.232O1619+01 

.24244122+01 

.61632565+01 

000041 

COVAR 

-97626021+01 

.36964213+01 

.35008504+01 

-.16000719-00 

.60677266+01 

000042 

COVAR 

.37219972+02 

.89507602-00 

.74552272-00 

.40220684-00 

.95687257-00 

000043 

COVAR 

.27100793+01 

-50692932+01 

.25373175+01 

.17656046+01 

-.46881713-01 

000044 

COVAR 

,31544700+01 

.14667230+02 

.13187752+02 



000049 

COVAR 

.95184916401 

.69113220+01 

.88646785+01 

.44309080+01 

.37431651+01 

OOOC46 

ec-VAR 

.3440Ci*r.7»2: 

.35:c7coo<oi 


2225*^4*^3 +*"** 

.31219773+01 

000047 

COVAR 

.74805143+01 

,70079423+01 

.42429308+01 

.39982355+01 

.95062911+01 

000048 

COVAR 

.48544640401 

.49499441+01 

.27374774+01 

.29444791+01 

.57372899+01 

0O0O49 

COVAR 

.50729561+01 

.29416253+01 

.29717698*01 

.19558660+01 

.19986400+01 

000050 

COVAR 

.34101902+01 

.31941114+01 

.32717162+01 

,56928901+01 

.61907175+01 

000051 

COVAR 

.35284390+01 

.30050691+01 

.69833051+01 

.57373951*01 

.39782751+01 

000092 

COVAR 

-95662936+01 

,43204475+01 

.50512070+01 

.27207938+01 

.33479249+01 

000093 

COVAR 

.53295736401 

.46932031+01 

.36034181*01 

.75018433+01 

.82316443+01 

000054 

COVAR 

,33267710+01 

.38098232+01 

.24236578+01 

.2651 7294+01 

.51118681+01 

000099 

COVAR 

.48265631+01 

.37793995+01 

.70766951+01 

.63806170+01 

.86976087+01 
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OODQ50 

COVAfl 26494707+01 - 

QOOQ57 

<OVAR 

.14704744-00 - 

000058 

COVAR 

.53009200*02 - 

000059 

COVAR - 

•,22509593*01 

000060 

COVAR -.20127919+01 

000061 

COVAR 

.91161002*01 

DOQ06Z 

COVAR 

.39052351*01 

000063 

COVAR 

.90973590*01 

000064 

COVAR 

.903051 72*01 

000065 

COVAR 

.15131915*02 

000066 

COVAR 

.81 560653*01 

000067 

COVAR 

.59410313*01 

000060 

COVAR 

.10960330*02 

000069 

COVAR 

.93032612*01 

000070 

COVAR 

.43050395*01 

000071 

COVAR 

.65711004*01 

000072 

COVAR 

.31010034*01 

OOOOT3 

COVAR 

.25031164*01 

ODOD74 

COVAR 

.53635000*02 

OD0075 

COVAR 

-.60693100*01 

000076 

COVAR 

-.11601605*02 

DOO077 

COVAR 

.70732130+01 

000070 

COVAR 

.45136420+01 

000079 

OWAR 

.89063820+01 

000080 

COVAR 

,03313652*01 

000081 

•WAR 

.14117970*02 

000002 

COVAK 

.93320693+01 

□00003 

COVAR 

.03905627*01 

CXXX304 

OOVAR 

.27528585*02 

000065 

COVAR 

,*7647301*02 

000086 

COVAR 

.15396030+02 

000007 

COVAR 

.25300704+02 

Dcxioea 

OOVAR 

.14244121+02 

000009 

COVAR 

.22493244+02 

000090 

OOVAR 

,95121075+02 

OOOD91 

COVAR 

.16004456+02 

000092 

OOVAR 

.34414124+02 

000093 

COVAR 

.43000513+01 

000094 

OCVAR 

.23645083+01 

030095 

COVAR 

.40020739+01 

000096 

OOVAR 

.40016750+01 

000097 

OOVAR 

.10031005+02 

000098 

COVAR 

.41551566+01 

000099 

COVAR 

.28342117+01 

DOOIOO 

COVAR 

,04690590*01 

OOOlOl 

COVAR 

.4*871692*01 

OOOIOZ 

COVAR 

.40576895*01 

000103 

COVAR 

, 9 1 6 / « 5 5 1 ^ i 

DD0104 

COVAR 

,49560077+01 

fXJOlOS 

COVAR 

.21281065+02 

000106 

COVAR 

.21452662+03 

000107 

COVAR 

,36359740*01 

000*08 

COVAR 

,16360267*02 

CXXI109 

COVAR 

,33940493*01 

OOOllO 

COVAR 

.30394506*01 

OOOtll 

COVAR 

,34706529*01 

fX>0tl2 

COVAR 

,17849460+01 


,76590752-00 

.29004777*01 

.14740525-00 

.23501945*02 

.66045053*01 

.40055606*01 

.94961002*01 

.91012062*01 

.60014969*01 

.02277669*01 

.67607252*01 

.57690475*01 

,62525052*01 

.65045769*01 

.61959697*01 

.20710610*01 

.76594085-00 

.10416247*01 

.25502506*01 

.26076915*02 

.65937751*01 

.43315008*01 

.99975557+01 

.95507591*01 

.50906703*01 

,95437159*01 

,90041533*01 

.10326594+02 

.17422523*02 

.17003708*02 

.10426264*02 

.14743792*02 

.16482672*02 

.10900059*02 

.17334939+02 

.34709601*02 

,27306034*01 

.19060092*01 

.40453460*01 

.43084301*01 

.26428291*01 

.32269907*01 

.31259096+01 

,26567734+01 

.40680440+01 

.42320746*01 

. 52553302*01 

.31544105*01 

.02790028*01 

.34650994*01 

,12624570*02 

.11183525*03 

.25994401*01 

.20761157*01 

,46314531+01 

.25869401+01 


-.16010528*01 - 

-.70749177*01 ' 

-.30246799*01 - 

-.07003701-00 ' 

.20961916*02 
.02326256*01 
.44924960*01 
.63655660*01 
♦54772540*01 
.56210070*01 
.62664110*01 
.12994063*02 
.65223772*01 
.59041790*01 
.45273624*01 
.16203404*02 
-.94750997-00 
-.15302163*02 
-.42929035*01 
-.34632219*01 
.24453757*02 
.069yJl263*OI 
.49940202*01 

.67359941+01 

,62005633+91 

.63134262+01 

.75400166+01 

.18665611+02 

.11044141+02 

.12573303+02 

.11564956+02 

.36561004+02 

.90199905+01 

,30606431+02 

.11674234+02 

,12055119+02 

.30140690*02 

.42106616+01 

.20096312+01 

.27010042+01 

.27385797+01 

.26290962+01 

,41364936+01 

.65373093+01 

.31176464+01 

.29523959+01 

.28606567+01 

.sr.o376Z6+01 

.26165750+01 

-.42543352+01 

,24480321+01 

.46594059+01 

.77303924+02 

.30617863+01 

.28601597+01 

.35623803+01 

.10047699+O1 


33224420+01 

,86234251*01 

,17496901+01 

,40709334+01 

,47323100+01 

.27606968+01 

.63051358*01 

,50472519*01 

.50423741*01 

.71519314*91 

.10407219*02 

.42061632*01 

.16019952*02 

.54767917*01 

.16604998*02 

.36734321*01 

,213C020Z*92 

.31550001*01 

,12824065*02 

.46457593*01 

.33001512*01 

.e950'.7961*0t 

.66909069*01 

.45411082*01 

.10912009+02 

.10224221+02 

.79355423+01 

.25320450*02 

.12408022*02 

,35929001+02 

.10013674*02 

,20524981*02 

.73607926*01 

.23043551*02 

.21120412*01 

.13551776+01 

,27045070*91 

,29434026*01 

,20345934*91 

.39005493*01 

.64847400+01 

.19793039+01 

.56737098+01 

,28217591*01 

.30612106*01 

-.46490452-01 

-.91409824*01 

.10001223*01 

-.17604962*01 

.21565959*01 

.20531092*01 

.37309828*01 

.21043453*01 


-.24777953*01 

-.30648Z36+01 

-.21301047+01 

-.45364016+01 

.41054335*01 
.39517795*01 
.16003301+02 
.13257359*02 
.39360126*01 
.06640605*01 
.77813252+01 
.50207100+01 
.19493933*02 
.92160462*01 
.18009220*02 
-.36075191*01 
-.15697665*02 
-.41 417114+01 
-.14913033+02 

.47161490+01 
.42900776+01 
.15971593+02 
.13493629+02 
.46625306*01 
.12441205*02 
.12909707*02 
.03170364+01 
.26206269*02 
.25537650*02 
.54530977*02 
.21706204+02 
.44914107+02 
♦70100087+01 
.22248609+02 

.10470369+01 
.21591266*01 
.76100952*01 
.63653137+01 
.20216123*01 
.41291670*01 
.41971929+01 
.24352960*01 
.57204728*01 
.67202036*01 
.11560847*02 
,65008956+01 
.20502022+02 
.s'l 73 1039-00 
-,4'^(743413*91 

.25345017*01 
.310110292+01 
.8004 *39 1+01 
.3945^158+01 



000113 

COVAR 

.42240608+01 

000114 

COVAR 

.56267141+01 

000115 

COVAR 

.44936745+01 

oootte 

COVAR 

.13406971+02 

oootir 

COVAR 

.77874967*01 

oooiie 

COVAR 

.20885105*01 

000119 

COVAR 

.34219839+01 

Dcxnzo 

COVAR 

-.13339377*02 

000121 

COVAR 

-.13251863+02 

000122 

COVAR 

.15790088+03 

000123 

COVAR 

-.15442987+02 

000124 

COVAR 

-.87529375+01 

000125 

COWAR 

.49165540+01 

000126 

COWAR 

,33452266+01 

000127 

COVAR 

.44058318+01 

00012B 

COVAR 

.50541981+01 

000129 

COVAR 

.56659295+01 

000130 

COVAR 

.36849940+01 

000131 

COWAR 

..46591196+01 

000132 

COWAR 

.14565491*02 

000133 

COVAR 

-86026637+01 

000134 

COVAR 

.49800909+01 

000135 

CCWAR 

.71898295+01 

000136 

COVAR 

-.86492268-00 

000137 

COWAR 

-.50989533-00 

000138 


.17327203+02 

000139 

COVAR 

-.14524A77+91 

000140 

COVAR 

-.16796910+01 

000141 

COVAR 

. 32550607-0* 

000142 

OOWAR 

.18139151+01 

000143 

CCWAR 

.11976308+01 

000144 

COVAR 

. 10077973+01 

0CO145 

COVAR 

,22610586+01 

000146 

COVAR 

,61779188-00 

000147 

OOWAR 

,67491376-00 

000148 

COVAR 

.33097997+01 

000149 

OOWAR 

.88500848-00 

000150 

COVAR 

.91205578-00 

000151 

COVAR 

.14620633+0* 

000152 

COVAR 

.11430746+01 

000153 

COVAR 

.12279202+01 

000154 

COVAR 

.77661507+01 

000155 

CCWAR 

.73212484-00 

□00156 

CC3VAR 

.10505424+01 

000157 

OOWAR 

.68706266+01 

000158 

COVAR 

,66030709*01 

OCW159 

COWAR 

.15630663+02 

000160 

COVAR 

,15201259*02 

000161 

COVAR 

.72338877+02 

000162 

COVAR 

.41799735:02 

000163 

COVAR 

.26728562+02 

000164 

COVAR 

.11890974+03 

000165 

COVAR 

.98330962*02 

000166 

COVAR 

.16157709*02 

oooter 

COVAR 

.86183519*02 

000168 

COVAR 

,14353229*02 

000169 

COVAR 

.81671889+02 



. 18502505+01 
.27769227 *^01 
.5223281 5+01 
.39994741+01 
.47508909*01 
,3t060 152*01 
.30405818+01 
-.18982635+02 
-.13568311+02 
-.01828662+01 
-.77197972+01 
.80553464+02 
,35723952+01 
,25022005+01 
.53944606+01 
.43250253+01 
.25037174+01 
.32807907+01 
.52472817+01 
.49408197+01 
.74982783+01 
.65610188+01 
.55528820+01 
-.13986376+01 
.96199417-01 
-.559CXI517-00 
-.49831312-00 
.39055138+01 
.12615095+01 
,62498125-00 
.12165587+01 
.11926076+01 
.78932373-00 
.77627830-00 
.56737077-00 
.46986055-00 
.11760384+01 
.62444388-00 
.10777113+01 
.82058982-00 
.13620929*01 
.16501724-00 
.66091383-00 
.40155540-00 
.80503303+01 
.63652312+01 
.27509669+02 
,29750730+02 
.10344833+02 

.34390049+02 

.13068625+02 

.65934391+02 

;33659419+02 

.63325221+02 

. 29527407+02 

.59687936+02 


.25434419+01 
.34712646+01 
.89684043+01 
.59548405+01 
,45006771+01 
.23975031+01 
.62507169+01 
-.12817145+02 
-.36761883+02 
-.11334143+02 
-.83564496+01 
.57769464+02 
.66244251+01 
.32616881+01 
.37548627+01 
.25601462+01 
,291 87314+01 
.36319976+01 
.88234575+01 
.68539010+01 
. 56424323+01 
.45143542+01 
.12795913+02 
-.50846461-00 
-.29078372+01 
-.94253881-00 
-.82782820-00 
.62561319+01 
.33219616+01 
.79204274-00 
.74486429-00 
.57662318-00 
.62739054-00 
.15472314+01 
.10544532+01 
.73557160-00 
.58983275-00 
.79787074-00 
.16277674+01 
.90170940-00 
.12468237+01 
.37032174-00 
.12413620*00 
.31083594+01 
.17179579+02 
.12201341+02 
.20399752+02 
.20878015+02 
.20139579+02 
.36709733+02 
, 77293060+02 
.27502782+02 
.48013902+02 
,26132204+02 
.12094319+03 
.22703419+02 
.10293521+03 


,21325010+01 
.45274682+01 
.53848447+01 
.41357053+01 
.11354019*02 
.25172587+01 
.55648507+01 
-.16108408+02 • 

-.35367139+02 
-.84552642*01 
-.21086223+02 

,27712584+01 

.21989867+01 

,39188325*01 

.31882508+01 

.23071127+01 

.51732802+01 

.73300628+01 

.43839127+01 

.13265269+02 

.51375642+01 

.13242022+02 

-.87849683-00 

-.34271111+01 

-.69481619-00 

-.23742799+01 

.11530214+01 

. 38161643-00 

.76317410-00 

.80004498-00 

.70775228-00 

.77417903-00 

.12339931+01 

.24606387-00 

.12109964+01 

.77740493-00 

.15384919+01 

.66720566-00 

.11263186+01 

.39511170-01 

.72520293-00 

.59722883+01 
.87287547+01 
.24139091+02 
.2708X744+02 
,14737105+02 
.19155549*02 
,87182449+02 
.20461577*02 
.90495088+02 
,34770442+02 
.99297449+02 
.29585771*02 
,79427103*02 


.21692167+01 
.64026065+01 
.48782720+01 
.49273030+01 
.11864442*02 
.46931487+01 
.56771081+01 
.25014025+02 
.14273030+02 
■.92627939+01 
•.20130738+02 ‘ 

.26550640+01 
.32165788+01 
.87499475+01 
.51250390+01 
.24442102+01 
.67187227+01 
.54515640+01 
.54075115+01 
.151261 74+QZ 
.84701067+01 
.15011698*02 
-.19069832*01 
-.1664-4834*01 
-.88467602-00 
-.27734318+01 

.69605637-00 
.13705122*01 
.31644600+01 
.13071675+01 
.48912475-00 
.86959226-00 
.72692125-00 
.65462425-00 
.21664174+01 
.12629289*01 
.36108294+01 
.73966274-00 
. 18567028+01 
,39132249-00 
.84082260-00 

.10002082+02 

-11986938+02 

.59919423+02 

.60707763*02 

.19183952*02 

.35931651+02 

.62261973+02 

.27469134+02 

.77933288+02 

,75321265+02 

,13741620+03 

.66133605*02 

.10954714+03 
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000170 

COVAR .11816456*03 

.91153413*01 

.19363242*02 

000171 

COVAR .44470342*02 

.3D277419+02 

.35527575+02 

000172 

COVAR .89153946*02 

.66374119*02 

.48951422*02 

000173 

CLA5SCS=S,C,0,W,R,A,Y, 

X,E 


000174 

PRIflT=lO MAX 01J=3Q000 

MD. SHCIW=03 


000175 

CMAM«:L5=01 ,09,11 ,12 



000176 

CMA7»JELS=06 , 09 , 1 0 , 1 1 



000177 

CHAI«CLS=01 ,06, 10,11 




. 14836825+02 

,65798096+02 


.18418329*02 

.30066092+02 



Input Deck Description 


Please refer to pages 61-64 as an example of a typical input deck. 


Card 1 


Field 1. n - the nuiriber of input channels (measurements). 
In this case, n=12. 

Field 2. m - The number of classes (populations). In this 
case, m=9. 

Field 3. k - The size of the "best" subset. The task is 
to pick the “best" k of n channels. In this 
case, k=4. 

Format* This card is read with the following format: 

(9X,I2,9X,I2,6X,I2) 


Cards 2-4 These cards contain the mean readings for the n channels of 
class 1. Field 1 is the mean response for the first channel, 
..., field 12 is the mean response for the twelfth channel. 
Three cards are required to hold the data for 12 channels. 
Format: (5X,5E15.8) 


Cards 5-28 These cards contain the mean readings for the n channels of 
classes 2 through m. In this case, three cards are required 

per class . 

Format: (5X,5E15.8) 


Cards 29-44 These cards contain the lower triangular covariance matrix 
for class 1. (Since the matrix is symmetric, it is only 
necessary to read the lower triangular portion). 


Field 1 - The variance of channel 1 for class 1 

Field 2 - The covariance of channels 2 and 1 for class 1 

Field 3 - The variance of channel 2 for class 1 

Field 4 - The covariance of channels 3 and 1 for class 1 

Field 5 - The covariance of channels 3 and 2 for class 1 


Field 


- The variance of channel n for class 1 


From this, it can be seen that there are 78 fields per class 
in the 12 channel case and that there are 16 cards/class. 
Format: (5X,5E15.8) 
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Cards 45-172 These cards contain the lower triangular covariance 
matrices for classes 2 through m. The fields are 
encoded in a similar to those for cards 29-44, There 
are 16 cards/class in this case. 

Format: (5X,5E15.8). 

Card 173 A symbol associated with each class. These symbold should 
be in the same order as the class statistics they represent, 
i.e,, symbol 1 represents class 1, symbol m represents 
class m. 

Format: (7X,36(1X,A1)). , 


Card 174 Field 1 . The number of top ranked channel sets to be 

printed for criteria, 1, 2, and 3, 

Field 2. The maximum pairwise divergence to be allowed 
by the program in its criteria tests. If a 
divergence is greater than this maximum, then 
It is set to it. 

Field 3. The number of specific channel combinations asked 
to be displayed, i.e., the number of “show" 
requests. If there are no show requests, this 
field should be set to zero. 


Format: (6X.I2,9X,I5,10X,I2) 

Cards 175-177 The specific channel combinations the investigator would 
like to have displayed. For example, the investigator 
has requested to see set (1, 9, 11, 12). There should be 
the same number of these cards as was specified in field 3 
of card 174. 

Format: (8X,24(1X,I2)) . 
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